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Abstract 

We study a critical behavior for the eigenvalue statistics in the two-matrix model in 
the quartic/quadratic case. For certain parameters, the eigenvalue distribution for one of 
the matrices has a limit that vanishes like a square root in the interior of the support. The 
main result of the paper is a new kernel that describes the local eigenvalue correlations 
near that critical point. The kernel is expressed in terms of a 4 x 4 Riemann-Hilbert 
problem related to the Hastings-McLeod solution of the Painleve II equation. We then 
compare the new kernel with two other critical phenomena that appeared in the literature 
before. First, we show that the critical kernel that appears in case of quadratic vanishing 
of the limiting eigenvalue distribution can be retrieved from the new kernel by means of a 
^ double scaling limit. Second, we briefly discuss the relation with the tacnode singularity 

in non-colliding Brownian motions that was recently analyzed. Although the limiting 
density in that model also vanishes like a square root at a certain interior point, the 
^ process at the local scale is different from the process that we obtain in the two-matrix 

(SJ model. 
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1 Introduction 

The two-matrix model in random matrix theory is the probability measure 

^ exp (-n Tr {V{Mi) + W{M2) - rMiMs)) dMi dMs, 

defined on couples (Afi, M2) of n x n Hermitian matrices. Here r > is the coupling constant, 
Zn is a normalizing constant and V and W are two polynomials such that the density is 
integrable. In this paper we are interested in the asymptotic behavior of the eigenvalues as 
n — cx). 

The two-matrix model is commonly assumed to be an effective model to generate new crit- 
ical phenomena for the eigenvalue statistics. In particular, it is believed to generate all (p, q) 
minimal conformal models, whereas in one-matrix models one only obtains the (p, 2) minimal 
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models |14 ^ 125 ^ [27]. However, even for non-critical bulk universality few rigorous results have 
been obtained and a general description of the asymptotic behavior of the spectrum remains 
an important open question. 

Recent progress has been made in \22\ [23l [Ml H2] where the authors considered the special 
case W{y) = + and analyzed the asymptotic behavior of the eigenvalues of Mi 

when averaged over M2. It turns out that the mean eigenvalue density for the matrix Mi has 
a limit as n — )■ 00, which can be characterized by a vector equilibrium problem. In particular, 
in [22] the authors showed that for the special case 

(1.1) W{y) = + ay^/2, and V{x) = ^x^, 
and the special values 

a = —1 and r = 1, 

the limiting eigenvalue density for Mi vanishes like a square root near the origin, which lies 
in the bulk of the spectrum, see also Figure [T] As a consequence, near the origin the local 
correlations are no longer in the universality class of the sine kernel and critical behavior is 
expected. It is the purpose of this paper to characterize that critical behavior. 

The eigenvalues of Mi, when averaged over M2, form a determinantal point process with 
a kernel that is expressed in terms of certain biorthogonal polynomials [281 133] • Our main 
result is a triple scaling limit of that kernel near the critical point. The limiting kernel is 
expressed in terms of the unique solution to a 4 x 4 model Riemann-Hilbert problem for 
which the associated ODE can be expressed in terms of the Hastings-McLeod solution to the 
Painleve II equation. This is the solution to the Painleve II equation 

(1.2) q"{a) = 2q{af +aq{a), 
characterized by 

(1.3) g(cj) = Ai(cr)(l + 0(1)), ascr-^+00. 

In the proof, we follow the approach of |24] and perform a Deift/Zhou steepest descent 
analysis for the Riemann-Hilbert problem characterizing the biorthogonal polynomials that 
integrate the two-matrix model. The crucial point in our proof is the construction of the 
local parametrix around the critical point in terms of the 4x4 Riemann-Hilbert problem 
that defines the kernel. To deal with the triple scaling limit, we will also need to modify the 
A-functions that are used in the normalization of the Riemann-Hilbert problem. 

The 4x4 Riemann-Hilbert problem that characterizes the critical kernel is an extension 
(we have an additional parameter) of a Riemann-Hilbert problem that appeared before in 
|20j . In that paper Delvaux, Kuijlaars, and Zhang analyze the local process near a tacnode in 
a model of non-colliding Brownian motions (for alternative approaches to the same problem 
see [Il|36j). In this situation, the limiting distribution of particles also vanishes like a square 
root in the interior of its support. However, we will show that at the local scale this process is 
different from the one we obtain for the two-matrix model. More precisely, the kernel that is 
of importance in the two-matrix model is constructed in a different way out of the extended 
(but for a certain choice of parameters same) Riemann-Hilbert problem. Note that this means 
in particular that the universality class is not determined by the vanishing exponent only. 

We recall that the Hastings-McLeod solution to the Painleve II equation also plays an 
important role in a different critical phenomenon in random matrix theory. It appears in 
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Figure 1: The limiting mean eigenvalue density for critical values a = —1, t = 1. 



the one-matrix model when the limiting density of eigenvalues vanishes quadratically in the 
interior of its support [U [12]. In that case the local behavior is expressed in terms of ^- 
functions that solve a 2 x 2 Riemann-Hilbert problem characterizing the Hastings-McLeod 
solution to the Painleve II equation. However, we stress that this 2x2 Riemann-Hilbert 
problem is essentially different from the 4x4 Riemann-Hilbert problem that we will use. 

It was already observed in [20j that the underlying Riemann-Hilbert problem for the kernel 
is intimately connected to the Hastings-McLeod solution of the Painleve II equation. In the 
context of the tacnode singularity in the Brownian motion model, this remarkable connection 
agrees with earlier results in |19) . In the critical case for the two-matrix model, we have an 
even stronger connection. Indeed, the ar- phase diagram tells us that it is possible to arrive at 
that situation, when we move away from the critical point along a certain curve. The second 
main result of the paper is that we can indeed retrieve the kernel associated with quadratic 
vanishing, by a double scaling limit of the new critical kernel. 
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2 Statement of results 

In this section we state our main results. The proofs will be given in the subsequent sections. 
We start by exploring some preliminary results on the two-matrix model that we will need. 
See also [24j and the references therein. 



2.1 Biorthogonal polynomials 

The two-matrix model can be integrated in terms of biorthogonal polynomials. These are 
defined as the two families of polynomials {pk,n}'kLo ^^'^ {lj,'n.}^oi where pk,n and qj^n are 
monic polynomials of degree k and j respectively, satisfying the orthogonality relation 

(2.1) / Pfc,„(x)g,-„(y)e-"(^(-)+^(^)---^) dxdy = 5,,hl^. 
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One can show that these polynomials exist and are unique |26j . Furthermore, they have 
properties that are typical for orthogonal polynomials: the zeros are real and simple [26] and 
they satisfy an interlacing property \22\ . The integrable structure of the polynomials also has 
been intensively investigated in [H Ul [5] . 

The correlation functions for the eigenvalues of Mi and M2 have a determinantal structure 
and can be expressed in terms of these polynomials |28[ H3] . Here we are only interested in 
the eigenvalues of Mi averaged over M2 and hence we content ourselves by describing the 
eigenvalues of Mi only. To this end, we define the transformed function Qj^n{x) as 



Q,, „(x) = e-'^^(^) / q,,n{y)e- 
JR 



-n{W{y)-Txy) 



dy. 



Then it is clear from (2.1) that pk.n and Qj,n satisfy the orthogonality condition 

Pk,n{^)Qj,ni^) = 0, j / k. 



The fact of the matter now is that the eigenvalues of Mi, when averaged over M2, form a 
determinantal point process with kernel 

n—l ^ 

(2.2) Kn{xi,X2) = -r^— Pk,n{^l)Qk,n{^2)- 

For more details on determinantal point processes we refer to [9} 134 1 [SS] , 140 1 1^1^ I46j . Note that 
in the language of [8], the eigenvalues of Mi form a biorthogonal ensemble. As the eigenvalue 
statistics are fully described by the kernel it remains to analyze the asymptotic behavior 
of this kernel near the critical point. 

Riemann-Hilbert techniques offer a natural approach to find the asymptotic behavior of 
Kn- We recall that the asymptotic behavior of orthogonal polynomials describing the one- 
matrix models (unitary ensembles), can be effectively computed using the Deift/Zhou steepest 
descent method |18j on an associated Riemann-Hilbert problem introduced in |30] . This pro- 
cedure was carried out in two important papers [16| ITT] in the regular situation. See also 
|15j for a detailed discussion. This development had a strong impact on the orthogonal poly- 
nomial literature. In many subsequent papers asymptotic results for (multiple) orthogonal 
polynomials were obtained through Riemann-Hilbert techniques. In particular, this method 
proved to be very successful to establish various universality results in random matrix theory. 

Also the biorthogonal polynomials describing the two-matrix model can be characterized 
by a Riemann-Hilbert problem 126 1 157 1 138) . Although the general situation is still open, for 
the special case of even V and W{y) = y^f-i + ay^/2, the steepest descent method has been 
successfully carried out for the Riemann-Hilbert problem [35] under the assumption that the 
limiting eigenvalue density is regular. See [23l I12j for the case a = and [24] for general a. 
We will proceed by discussing some of these results. 

2.2 Limiting mean density 

Prom this point we assume that V{x) is an even polynomial with positive leading coefficient 

T^(y) = ^+a|-. 



5 



For technical reasons we will also assume that n = mod 3. It then follows from the results 
in [23] that the mean eigenvalue density for Mi has a limit as n — )• oo, i.e., there exists an 
absolutely continuous measure /ii on M such that 

(2.3) liin -Knix,x) = ^. 

n-5>oo n ax 

Moreover, from the results in [2l] one also deduces that fii is the limiting normalized counting 
measure on the zeros of pn n 



n ^ — ' 

X ■■ Pn,n{x)=0 

This was also established in [22] for the special case V{x) = x^/2 using different methods. A 
crucial point in the analysis of [24J is that the limiting measure //i can be characterized as 
the first component of a vector of measures (//i,/i2,^3) that is the unique minimizer of an 
energy functional E, which we will now discuss. We need some basic notions from potential 
theory, for which we refer to the standard work [45J. 

Let us first introduce the logarithmic energy for a measure /i on C 

HfJ') = / / log - — - — r dfi{x) d/i(y). 
JJ \x-y\ 

For two measures fi and i' the mutual logarithmic energy is defined as 

Iif^,^)= II log] — - — rd/i(x)di/(y). 
J J \x-y\ 

Then the equilibrium problem is to minimize the energy functional 

3 2 

(2.5) E(i/i, 1^2,1^3) = ^lii^j) - + / Viix)di^iix) + I V^{x)dv^{x), 

3=1 i=i •' 

among all measures ui, U2, and that satisfy 

(a) the measures have finite logarithmic energy; 

(b) ui is a measure on M with z^i(M) = 1; 

(c) U2 is a measure on iM with V2{iM^ = 2/3; 

(d) is a measure on M with z^3(M) = 1/3; 

(e) f2 < 172 where a2 is a certain measure on the imaginary axis. 

The vector equilibrium problem clearly depends on the input data Vi, V3, and (T2 that we 
describe next. We will only deal with the situation a < 0, which is the situation relevant to 
us. For a > we refer to [24J. The external field Vi that acts on vi is given by 

Viix) = Vix) + min(Ty(s) — txs). 
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Now define 



(2.6) x\a) = -' 

T 

Then, the function M — )• M : s i— )• W{s) — txs has a local maximum and a local minimum 
besides the global minimum if and only if x G (— x*(a), x*{a)). In that case the external field 
V3 is defined as the positive difference between these local extrema. For all other real values 
of X we define V^{x) = 0. Finally, to describe the measure (T2 that acts as a constraint on V2 
we consider the equation 

+ as = Tz^ with z G iM, 
and let s{z) denote its solution with the largest real part for z G iR. Then 

d<J2(2;) r 



\dz\ TT 



Re s{z). 



and the support of o"2 is zM. 

In [24j it is proved that the above vector equilibrium problem has a unique minimizer 



(/ii, /i2, /^s) (see also [33]) and that (2.3) holds. It is also proved that each has an analytic 



density with respect to the Lebesgue measure. 
2.3 Phase diagram 

From now on we also assume that V{x) = ix^, so that we are in the situation given by ( 1.1 ). 



In order to get a better understanding of the nature of the critical point that we will consider 
in the present paper, we will first discuss the behavior of the supports of the measures fii, 
— fJ'2, and and how they depend on a and r. As proved in [22l |2l|, the supports of the 
measures /ii, (J2 — ^2^ and ^3 have the following form 

supp(Aii) = [-/3o,-/3i] U [/3i,/3o], 
supp(fT2 - ^12) = \ [-1^2, 1(^2), 
supp(//3) =M\(-/33,/33), 

for some /3o > /3i > 0, /32, /33 > that all depend on the values of a G M and r > 0. We 
distinguish a number of cases, depending on whether /32, or /33 are equal to zero, or not. 
At least one of these is zero, and generically, no two consecutive ones are zero. The situation 
in summarized in the phase diagram in the ar-plane shown in Figure [2| 

Case I: /3i = 0, /32 > 0, and = 0. Thus, in this case there are no gaps in the supports of 
the measures ^1 and /i3 on the real line. The constraint (J2 is active along an interval 
[— z/32, on the imaginary axis. 

Case II: /3i > 0, /32 > 0, and (3^ = 0. In Case II there is a gap in the support of /ii, but 
there is no gap in the support of /i3, which is again the full real line. The constraint is 
active along an interval along the imaginary axis. 

Case III: /3i > 0, /32 = 0, and /33 > 0. In Case III there is a gap in the supports of /xi and 
/X3, but the constraint on the imaginary axis is not active. 
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Case IV: /3i = 0, /32 > 0, and /^s > 0. In this case the measure fxi is stih supported on 
one interval. However there is a gap (— /33,/33) in the support of /is. As in Case I, the 
constraint (T2 is active along an interval [— i/32,i/32] on the imaginary axis. 

In Figure [2] we plotted a phase diagram that shows which values of (a,T) correspond to the 
different cases. The different cases are separated by the curves given by the equations 



r = Va + 2, — 2 < a < 00, and r = \ , —00 < a < 0. 

V a 

On these curves two of the numbers /3i,/32, and (3^ are equal to zero. For example, on the 
curve between Case III and Case IV, we have /3i = /32 = 0, while (3-^ > 0. Finally, note the 
multi-critical point (a, r) = (—1, 1) in the phase diagram, where /3i = /32 = /^s = 0. All four 
cases come together at this point in the ar- plane. This point has our main interest. 

As long as we consider points (a, r) that are not on the curves, the local correlations are 
governed by the sine kernel in the bulk of the spectrum and the Airy kernel at the edge of 
the spectrum. Critical phenomena occur at the curves that separate the different cases. 

When we cross the line r = \/a + 2 for a > —1, the support of //i turns from two intervals 
into one interval. On the critical curve, the intervals meet at the origin and there the density 
vanishes quadratically. This indicates that in a double scaling limit the kernel converges to a 
kernel that is constructed out of ^'-functions for the Hastings-McLeod solution of the Painleve 
II equation [6l[l2l[l3]. 

When crossing the line r = sj —Xja for a < — 1, we again have a transition of two intervals 
merging at the origin. However, due to the fact that a — [12 also closes simultaneously, the 
vanishing at the origin occurs with an exponent 1/3. This indicates that the local correlations 
are governed by the Pearcey kernel (TJ HOl [HI SH HZ] • 

The other transitions, represented by the dashed lines in Figure [2l do not concern ix\. 
They take place on the non-physical sheets of the spectral curve and, therefore, they do not 
influence the local correlations of the eigenvalues of M\ , which are again described by the sine 
and Airy kernels. 

In the next paragraph, we will describe the limiting process near the point (a, r) = (—1,1), 
where there is a simultaneous transition in the supports of all three measures /ii, cr — /i2 and 

fJ'3- 

2.4 Limiting kernel 

We now come to the main results of the paper. We study the process near the critical 
parameters r = 1 and a = —1, by means of a triple scaling limit. To this end, we rescale a 
and r near the critical values in the following way 

,2,7, (^) (J) 

for a, 6 G R. Note that the vectors (2 l) and (—1 2) are respectively tangent and normal 
to both critical curves in the point (a, r) = (—1, 1). We also scale the space variables with 

X = un~'^^^, and y = vn~'^^^, 

and compute the limiting behavior of Knix, y) as n — )• 00. 

The limiting kernel is characterized by the solution to a Riemann-Hilbert problem (= RH 
problem) that we will first discuss. 
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Figure 2: The phase diagram in the ar-plane: the critical curves r = ^Ja + 2 and r = y — ^ 
separate the four cases. The cases are distinguished by the fact whether is in the support 
of the measures Hi, 02 — lJi2, and )U3, or not. 

RH problem 2.1. Fix parameters ri, r2, si, S2; and t. We search for a 4 x 4 matrix-valued 
function M(Q satisfying 

(1) M is analytic for C G C \ Em; 



(2) M+(C) = M_(C) Jfe, forCeTk,k = 0,..., 9; 



(3) As C ^ 00 with C G C \ Sm t^e have 



(2.8) M(C)=(/ + 0(r'))5(C)^ 



X 



diag ( e-^^«)+tC^ e-Vi(C)-tC^ eV2(C)+tC^ gVUO-tC 



where 



(2.10) 



(2.9) 



V'i(C) 



V^2(C) 



^nC^/^ + 2.iC^/2, 

lr2i-Cf/' + 2s2{-C)'/' 



and 



(2.11) 



x/2 



1 



/I -f 0\ 

1 i 

-i 1 

\0 i ij 



and, finally, 



(2.12) 




(4) M{Q is bounded near C = 0. 
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Figure 3: The jump contour Sjv/ in the complex (^-plane and the constant jump matrices 
on each of the rays r^, A; = 0, . . . , 9. 

The fractional powers in are chosen such that these 

maps are analytic in C\ (— oo,0] and taking positive values on the positive part of the real 
line. The fractional powers in Q ^ (— C)'^^^> C ^ (~C)^^^ and Q i— )• (— C)^^'''' are chosen such 
that these maps are analytic in C\ [0, oo) and taking positive values on the negative part of 
the real line. 

The contour T^m is shown in Figure\^ and consists of 10 rays emanating from the origin. 
The function M{(^) makes constant jumps on each of the rays Tj^. These rays are deter- 
mined by two angles (pi and (p2 satisfying < 921 < 992 < 7r/2. The half-lines F^, /c = 0, . . . , 9, 
are defined by 



and 



Fi 

- 



„i(7r-(p2)m + 



F2 



k = 0, 



All rays are oriented towards infinity. 



Notations hke M±{z) are used throughout the paper and should be understood as follows. 
Suppose that z belongs to an oriented contour (here z E Sa/)- Locally the contour divides the 
complex plane into two parts. The part on the left (right) side when traversing the contour 
according to its orientation will be called the + (-) side. Then M+(z) (M_(z)) denotes the 
limiting value of M{z) as z approaches z from the + (-) side of the contour. 

In case t = 0, the RH problem 2.1 has already appeared before in connection with the 



tacnode singularity in the model of non-colliding Brownian motions [20] . The following result 
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on the solvability of the RH problem for M was already obtained in |20j for the special 
situation t = 0. Here we will also need the statement for t S M. 



Theorem 2.2. Assume ri = r2 > 0, si = S2 ^ M, and t € M. Then RH problem 2.1 has a 
unique solution. 

The proof of this theorem will be given in Section [5] We remark that we do not prove 
the existence by a vanishing lemma, which does not seem to work for i £ M. Instead, we 
prove the existence by combining the existence for t = 0, as obtained in [20], with standard 
considerations from the theory of isomonodromy deformations. 

As mentioned in the introduction RH problem 2.1 is connected to the Hastings-McLeod 
solution to the Painleve II equation, see (1.2) and (1.3). This will be made more precise in 
Section 5.1 when we discuss the Lax pair behind the RH problem. For now, we only mention 
that it is possible to retrieve the Hastings-McLeod solution by considering the following limit 



= i2-'l\ 



lim ( CM(C) diag ( e'^^CO-tC^ ^^liO+K ^ ^-MO 

"\[ ~/\2s-t'' 

where q stands for the Hastings-McLeod solution and we assumed ri = r2 = 1 and si = S2 = s. 

Now that we have established the existence of the solution of the RH problem, we will 
define the limiting kernel. It is convenient to first transform M in the following way 

M{C) = MiiC)-''C±, 



(2.13) 
where 

(2.14) 



iblmC > 0, 



C4 



/I 













-1 

















-1 


\o 





-1 


0/ 



and C_ 





-1 





^\ 


-1 

















-1 





Vo 








1/ 



Naturally, the function M is also characterized by a RH problem, that can be easily obtained 
out of the RH problem for M. The precise formulation of that RH problem will be given in 
Section 4.8, where it is used in the construction of the local parametrix. 
We now define Kcr by 

(2.15) Kcr{u,v;s,t) 



1 



27ri(n 



(-1 1 0) M~^{v)M+{u) 



n\ 

1 





where M_|_ is the -|— boundary value of M, as in (2.13), with M(u) the unique solution M to 

1 > 0, 



RH problem 2.1 with parameters 
(2.16) 



ri = r2 

Sl = S2 

t e M. 



s G 



A straightforward calculation shows that kernel K{u, v) depends analytically on u, v, also in 
the origin. 

The following theorem is the main result of this paper. 
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Theorem 2.3. Let he the kernel in (2.2) describing the eigenvalues of Mi when averaged 



over M2 with V and W as in (|1.1). Set 

-1' 



+ an 



-1/3 



1 y ' \i 

for a, 6 G M. Then for n — ?• 00 and n = mod 6 

1 



+ 



lim 



n— >oo n 



2/3 



Kn 



uniformly for u,v in compact subsets o/M, where Kcr is as in (2.15) and (2.16). 



To the best of our knowledge, the kernel Kcr and the associated process have not appeared 
in the literature before. 



Remark 2.4. The assumption n = mod 6 in Theorem 2.3 is mainly for technical reasons. 
We recall that also in \23\ [2^ H2] the assumption n = mod 3 was used. Here we will need 
the additional assumption that n is even. However, we believe that these assumptions are for 
technical reasons only and that with some additional argument this condition can be dropped 
entirely. 

2.5 A double scaling limit for 

The nature of the critical point suggests that it is possible to take interesting limits of the 
kernel. Indeed, if we walk away to the right starting from the critical point and staying on the 
curve r = \/a + 2, we expect to retrieve the Painleve II kernel that describes the transition 
from Case I to Case II. In the local variables a, b this corresponds to setting b = — and 
taking the limit a — )• +00. This leads to t = — a and s = as parameters for Kcr- In 
fact, to obtain the limit in the most general form one needs to perform a double scaling limit. 
Hence we set 



1„2 



(2.17) 



(l-o^)' 



for fixed o" € M and take the limit a — )• +00. 
Consider the following RH problem. 

RH problem 2.5. We look for a 2 x 2 matrix-valued function ^(C) satisfying 

(1) ^'(C) is analytic for ( £ C \ T,^; 

(2) +{C) = ^-iC) Jk, for C e Tk, k = 1, ... ,4; 

(3) As ( ^ 00 we have 



^{C) = {i + o{c')) 











(4) ^(C) is bounded near C = 0. 
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Figure 4: The jump contour T,^ in the complex (^-plane and the constant jump matrices 
on each of the rays r^, A; = 1, . . . , 4. 



The contour Sij, = Fi U r2 U U where 



2- 



ri = e'^*/^M+, r2 = e^^'^/6M+, = -Fi, r4 = 
All rays are oriented towards infinity. See also Figure^ 

This RH problem was introduced by Flaschka and Newell in pS] . They showed that from 
this RH problem one can retrieve the Hastings-McLeod solution for the Painleve II equation. 
More precisely, define q{y) by 



(2.18) 



q{v) = lim C*i2(C;i^)e' 

s>oo 



where ^'12 is the 12-entry of ^, then q solves the Painleve II equation (1.2) and satisfies (1.3). 
One can show that there exists a unique solution ^ to RH problem 2.5 if and only if the 
Hastings-McLeod solution q has no pole at v. Since it is known that this solution has no 
poles on the real axis [32], it follows that ^ exists for all G R. 
We now define the kernel K-pn by 



(2.19) 



1 



27ri{u — v) 



[1 -l)^-\u;,^)^{v;u)(^^ , 



where ^'(C,z^) is the unique solution to RH problem 2.5 The kernel Kpn is the kernel that Are u,v 

appears as universal object in random matrix theory, when the limiting mean density vanishes inter- 

quadratically at an interior point of the support, see [6l[l2]. changed? 
Then we have the following limiting behavior of the critical kernel. 



Theorem 2.6. Let Kcr be as defined in (2.15) and Kpu as in (2.19). There exists a function 
h such that 



h{y, a) 



lim 2^/^a :r'7 Kcr ( 2^1^ ax, 2^'^ ay; y , -« ( 1 



a 



i^Pll(x,y;25/V), 



uniformly for x, y in compact sets. 



2.6 Comparison with tacnode situation in non-colliding Brownian motions 

As mentioned in the introduction, the phenomenon of square root vanishing of a limiting 
density at an interior point of its support can also be obtained in a model of non-colliding 
Brownian motions [H [201 [36] . 



13 



Consider 2n particles on the real line, half of which start at time S = —1 in position 
a and the other half start in position —a. The walkers perform Brownian motions and are 
conditioned to return to their starting positions at time 5 = 1. Moreover, we condition the 
particles never to collide. For re — )■ oo and a proper scaling of a, the paths of the particles 
in the space-time diagram fill out two ellipses that touch at 5 = 0, thus forming a tacnode. 
The limiting density of particles at = is described by two touching semi-circles. Hence, 
we see square root vanishing at the origin, precisely as in the multi-critical situation in the 
two-matrix model that we just discussed. In |20j the authors show (in a more general setup) 
that in a double scaling limit the process at the local scale near the tacnode is governed by 
the kernel 

n\ 



(2.20) 



Kti,c{u,v;ri,r2,si,S2) 



1 



2Tri{u — v) 



(-1 1 0) M-^{v)M+{u) 




1 



for u,v > 0, where M is the solution of the RH problem 2.1 Here rj and sj are expressed in 
terms of the relevant parameters of the problem and t = 0. The kernel for general w G M is 
defined by analytic continuation. 

Although the critical situations in the two different models share the same vanishing 
exponent, the kernels -fCtac and Kcr differ. This may not be immediately obvious from the 



representations (2.15) and (2.20), but can for example be seen by computing the asymptotic 



behavior along their diagonals. Based on RH problem 2.1 one shows that this asymptotic 
behavior is different for both kernels, which proves that the processes are indeed different. 
Indeed, for the kernel in the tacnode situation we have 

S2 1 



(2.21) 



Ktac{u,u;ri,r2,si,S2) 



r2yJU 



TT 



7^^/^l 



Attu 



Ree2'^2H + 0(n-3/2) 



as n — 7- +c«, where we recall the definition of V'2 (2.10). Note that there is no constant term 



in the expansion and that 1/u has a highly oscillatory coefficient of modulus 1 . For the kernel 
appearing in the multi-critical situation in the two-matrix model we have 



(2.22) 



Kcr{u,u;s,t) 



'u t s 
+ - + —=- 



as ti — )• +00. Here we do have a constant contribution, but there is no 1/u term. Therefore, 
even for well-chosen parameters, these two expansions cannot be the same and hence the 
kernels are essentially different. 



The computations that lead to (2.21) and (2.22) follow by rather straightforward com- 



putations on the asymptotic behavior of the function M and are given in Appendix A for 
completeness. 



2.7 Overview of the rest of the paper 

The rest of the paper is devoted to the proof of our main results. The proof of Theorem 



2.3 is based on a Deift/Zhou steepest descent analysis on a RH problem characterizing the 
biorthogonal polynomials. This RH problem will be introduced in Section 4. The steepest 
descent analysis will be given in the same section. In this analysis we will make use of certain 
meromorphic functions on a Riemann surface that we will define first in Section 3. In Section 



5 we prove Theorem 2.2 The proof of Theorem 2.6 can be found in Section 6 



14 



3 Meromorphic functions and Riemann surfaces 



In this section we define several meromorphic functions and derive some of their properties 
that we will need later on. In particular, we define the A-functions that we use in the nor- 
malization step in the steepest descent analysis. The definitions as presented here differ at 
some point from [2l], except in the critical situation a = —1 and r = 1. To emphasize this 
difference we speak about modified Riemann surface and modified ^, A-functions. These mod- 
ifications are typical when dealing with double scaling limits, see for example [71 
for the case of more than two sheets. 

Throughout this section, we will introduce new functions {9, w, ^, A, h, g, H, G, . . .) 
and constants (7, c, . . .). Although it is not indicated in the notation, these functions and 
numbers depend on a and r. We will always assume that the parameters r and a are close to 
their critical values r = 1 and a = —1. When dealing with functions or numbers associated 
with the critical parameter values we add a star to the notation, thus, we write w* , .^*, A*, 
7 , c ,. . .. 

3.1 Definition of the functions 6j 

We start with defining a function ^ on a three-sheeted Riemann surface. The definition is 
exactly the same as in |24j . 
Consider the equation 

+ as = TZ, 

for a < 0. If 2; = xEM and \x\ < x*{a) = ^ then the equation has three real 

solutions Sj = Sj{x) that we take to be ordered such that 

W{si) — TXSl < W{S2) — TXS2 < W{s3) — TXS3. 

In other words, the function W{s) — txz has a global minimum at si, a local minimum at S2 
and a local maximum at S3. The functions x 1— )• Sj{x) extend to a meromorphic function in 
the following way. 

Consider the three-sheeted Riemann surface S with sheets 

(Si =C\iR 

(3.1) < S2 =C\((-oo,-x*(a)]U [x*(a),oo)UiM) 
i cSs =C\((-oo,-x*(a)]U[x*(a),oo)) 

Then the function sj has an analytic continuation to Sj, which we also denote Sj. It is 
straightforward to check that S2 and S3 can also be obtained by analytic continuation of si 
onto ^2 and ^3. 

Finally, we define the functions 6j by 

(3.2) ej{z) = -W{sjiz)) + Tzsj{z), j = 1,2,3. 
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We will need the asymptotic behavior as z — )• oo that is given by |24^ Lem. 2.4] 



(3.3) ei{z) 





^(tz)'^/^ + ^ - 


54 


^ ; T \^ J , 




in 


lU IV, 




- f a;2(Tz)2/3 + 


a2 
6 


- fJc^(rz)-2/3 + O {z- 


4/3)^ 


in 


II, 




- f L^(tz)2/3 + 


a2 
6 


-fJc.2(rz)-2/3 + o(z 


-4/3) 


in 


III, 


[HrzY"- 


- f a;2(Tz)2/3 + 


a2 
6 


- fJt^(rz)-2/3 + O (z- 


4/3)^ 


in 


I, 


< lirzf/^- 




54 


;rz)-2/3 + O , 




in 


II U ///, 




- f w(tz)2/3 + 


6 


-f^^2(^^)-2/3 + 0(^ 


-4/3) ^ 


in 


IV, 


{ |l^2(^^)4/3 


- f a;(r2)2/3 + 


a2 
6 


^oj^{tz)-^/-^ + 0{z 


-4/3) ^ 


in 


lull, 


\Hrzf'- 


- f a;2(rz)2/3 + 


a2 
6 


^u(rz)-y^ + 0{z- 


4/3) ^ 


in 


III U IV. 



(3.4) e2iz) 



(3.5) e^iz) 



Here, I, //, ///, and IV denote the four open quadrants of the complex plane and 



3.2 Modified Riemann surface TZ 

Let us first introduce an auxiliary parameter 7 that is completely determined by a and r but 
will prove to be convenient for notation. We define 7 = 7(0, r) as the solution of 



ar 



2/3 _ ± 
7 



97' + 5r4/3 



7, 



(3.6) 

that tends to 1 as r — )■ 1 and a — )• — 1. In the triple scaling limit, i.e. we let a and r depend 
on n as in (2.7) while n — )• 00, we have 



1 + |an 



(3.7) 7 

Next we introduce a four-sheeted Riemann surface TZ with sheets 

7^l = c\ [-c,c], 7^2 
7^3 = c\(Mu«), 7^4 



\([-c,c]UiM), 
\M, 



where 
(3.8) 



16 ^ 
3^/3 



3/2 



We connect the sheets TZj to each other in the usual crosswise manner along the cuts [— c, c], 
M, and zM. For example, T^i is connected to 7^2 along the cut [— c, c] and so on. Moreover, 
the Riemann surface is compactified by adding two points at infinity. The first point ooi is 
added to the first sheet, while the second point 002 connects the other sheets. This Riemann 
surface TZ has genus zero. It is shown in Figure [5| 
Let us now consider the equation 



zw 



(3.9) {w^ + j^'f 

We claim that this equation defines an algebraic function on the Riemann surface TZ. Indeed, 



for each z G C equation (3.9) has four solutions, denoted Wj{z), j = 1,2,3,4. We order the 
solutions such that 



(3.10) 



\'Wl{z)\ > \w2iz)\ > \W3{Z)\ > 1^4(2)1. 
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Figure 5: Modified Riemann surface for values of the parameters a and r close to a = — 1 
and r = 1. Note the points at infinity: ooi belongs to TZi, while 002 is common to the other 
sheets. 



For certain z € C this can be done in different ways. In that case we arbitrarily pick a 
particular order. It follows from the next lemma that the functions Wj{z) can be interpreted 
as a meromorphic function on the Riemann surface IZ. 

Lemma 3.1. The function Wj{z) is analytic on TZj, j = 1, 2, 3, 4, and the following relations 
hold 







[x] 


= Wj-{x), 


X e (-c,0)U(0,c), j = l,2, 






[x] 


= W2,^{x), 


X £ (-c,0) U (0,c), 






{< 


= 


zGiM\{0}, j = 2,3, 




'W2,± 


i< 




z G iM\{0}, 






[x] 


= Wj-{x), 


xGM\{0}, j = 3,4, 




W3,± 


[x] 


= W4^^{x), 


X G M\ {0}. 


Here [— c, c] 


and 


R 


are oriented from left to right and iM. from —ioo to ioo 


function 




— > 


C : TZj 9 z 1— )• Wj{z) has an analytic continuation to a 


function w : 




C. 


Moreover, we 


have the symmetry condition 


(3.11) 






Wj{z) = - 


-Wj{-z), zGC\(MUiM). 



Proof. For the case 7 = 1 the relations on the cuts follow from [22, Th. 5.3] and the proof of 
that theorem. A rescaling proves the relations for general 7. 

Due to symmetry both {wj{z),j = 1,2,3,4} and {—Wj{—z),j = 1,2,3,4} solve the al- 
gebraic equation (3.9). The fact that both sets of solutions have to satisfy (3.10) proves 
( [3lT| ). □ 

The following lemmas on the asymptotic behavior of Wj{z) for z — ?• and 2; — )■ cx) will be 
important to us. Recall the notation uj = e^'^*/^. 
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Lemma 3.2. The meromorphic function w on TZ has the following asymptotic behavior as 







W2{z) 



wsiz) 



nj,{z) = ir-l'' + le-'lW'^z^l'' + iz + ^^e 



'i^m + le-3.V4^3/4^1/2 ^ 1^ ^ ^g3.V4^-3/4^3/2 ^ ^(^2) 
i^3/2 + lg,rV4y/4^1/2 + 1^ ^ _^g-.*/4^-3/4^3/2 ^ ^(^2)^ 

'i^3/2 ^ lgW4^3/4^1/2 ^ 1^ ^ _^g-W4^-3/4^3/2 ^ ^(^2)^ 
j^i^3/2 + lg-3«/4^3/4^1/2 ^ 1^ ^ ^^3.^/4^-3/4^3/2 + ^(^2) 

= i^3/2 ^ 1^-3.^/4^3/4^1/2 ^ 1^ ^ ^e3../4^-3/4^3/2 ^ ^(^2)^ 



m / U //, 

in III, 
in IV, 

in III, 
in IV, 

in lU II. 



In the remaining quadrants the expansion is determined by (3.11) and from above. 

Lemma 3.3. The meromorphic function w on TZ has the following asymptotic behavior as 
z — )■ oo 



wi{z) = z — 2"f' 

>,2 



W2iz) 



wsiz) 



W4^{z) 



31 



+ Oiz- 



m 



TZ + 173^-1 + ^7^2-5/3 + 0{z-y^), inlU IV, 

^2^2^-1/3 + 2^3^-1 + 7^^4^-5/3 + 0(^-7/3)^ jj^ 



ajj z 



1/3 + 2^3^-1 + 7^2^^-5/3 + 0(^-7/3)^ jjj^ 



r ^2^2^-1/3 + 2^3^-1 + 7^^4^-5/3 + 0(^-7/3)^ 



1 + 7^4^-5/3 + 0(^-7/3)^ 



in // U ///, 

in IV, 



^2^-1/3 ^ 173^;- 

(^^^2^-1/3 + 2^3^-1 + 7^2^^-5/3 + 0(^-7/3)^ 

[,^^2^-1/3 + 2^3^-1 + 7^2^4^-5/3 + 0(^-7/3)^ j y jj^ 

1 ^2^2^-1/3 + 2^3^-1 + 7^^4^-5/3 + 0(^-7/3)^ jjj y IF. 



3.3 Modified ^-functions 

We define the modified ^-functions for j = 1, 2, 3, 4 as 

Wj{z)^ + (373 - l)Wj{z)^ + T^/'^J^ 



(3.12) 



Wj{z){wj{z)^ + 73) 
Then, the analytic extension ^ : 7^ — t- C of the function 



z E TZj. 



(3.13) 



■.TZjBzt-^ ij{z), 



to the full Riemann surface TZ is meromorphic on 7^. This is the content of the following 
lemma. 



Lemma 3.4. The function ij{z) is analytic on TZj, j = 1, 2, 3, 4, and the following relations 
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hold 



(3.14) 




z G (-c,0) U (0,c), 


j 


(3.15) 


6,±(^) = 6,t(^)' 


z G (-c,0) U (0,c), 






6,+ (^) = -^2,-{z), 


zGi]R\{0}, i 


= 2,3 






z G iM\ {0}, 






= 6-(^), 


zGM\{0}, j-- 


= 3,4, 






z G M\{0}. 




Moreover, 


we /ia?;e the symmetry condition 






(3.16) 




z G C \ (RUiM). 





Proof. This is immediate from Lemma 3.1 and (3.12). 



□ 



Combining (3.12) with Lemma 3.2 and Lemma 3.3 leads to the asymptotic behavior of 
the ^-functions. 

Lemma 3.5. For z in a neighborhood of zero, we have 
Uz) = Cz-^/^+giz) + z^/^hiz) 



(3.17) 



in lU II, 

-Cz-^/^ + g{z) - z^/^h{z) in III, 
Cz-^l"^ + g[z) + z^l'^h{z) in IV, 

Cz-^l'^ + g{z) + z^l'^h{z) in III, 
-Cz-^l^ + g{z) - z^l^h{z) in IV, 



-Cz-'^'^ + g{z)- z^/^h{z) 



in lU II, 



where 
(3.18) 



(J ^ g37ri/4^1/4 



and g and h are analytic functions in a neighborhood of zero with 



X _ 54/3 
47 4 ' ' 



m 



TTi/i -^/A ( 3 2 



7 



3^2 _ ^ , 3 4/3. 
4 ' 87 ^ 8 ' 



In the remaining quadrants the behavior is determined by (3.16) and (3.17). 
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Lemma 3.6. The meromorphic ^-function has the following asymptotic behavior as z ^ oo 



1 -^4/3^5 _ 3 

- + f- 

Z Z'^ 



+ 0{z-^), 



in 



r ^4/3^1/3 _ I + I -I + 0(z-5/3), 



+ 0(z-5/3), 



^4/3^^1/3 . 


- iaT^/3^2^ 


-1/3 


+ 


3^ 


^4/3^2^1/3 


— |ar^/3(^2 


-1/3 


+ 


^Z 
3^ 


^4/3^^1/3 . 




-1/3 


+ 


^Z 
3^ 


T-4/3^1/3 _ 


1^^2/3^-1/3 


+ \ 


Z~ 




^4/3^2^1/3 


— ^ut'^/^ujz 


-1/3 


+ 


^Z 
3^ 


^4/3^2^1/3 


- ^ar'^/^ojz' 


-1/3 


+ 


^Z 

3^ 


^4/3^^1/3 _ 


- \aT^'^oj^z 


-1/3 


+ 


^Z 

3^ 



+ 0{z-^l^), 



in / U IV, 
in II, 
in III, 

in I, 

in II U ///, 
in IV, 

in lU II, 
in III U IV. 



Remark 3.7. Note that this behavior coincides with that of the unmodified ^-function in [231 
Lem. 4.f0] up to and including the 1/z term. ActuaUy this motivates the definitions of ^ 
( [3I2I ) and 7 ( [3^ . 

As a corollary to both lemmas we see that has simple poles in both points at infinity. 
^ also has at most simple poles in the two points at zero. Only if C = C(a,r) = 0, the 
^-functions are bounded in zero. Note that in the multi-critical point a = — 1, r = 1 the 
function ^ has simple zeros in the two points in zero. 

We will also need the following integrals of ^j. 

Lemma 3.8. We have the following definite integrals 



(3.19) 
(3.20) 



Im J ^i-t (x) dx = ±7r, 
Im J S,2±ix)dx = =F7r, 



with c as in (3 A 



Proof. In view of (|3.15|) the expressions (|3.19|) and (|3.20|) are equivalent. Hence, we only 

1 



prove the first one. Note that by (3.14) 



Im / £,i,±{x) dx 



1 

"2i 



(6,+ -6,-)(x)dx 



2i 



^i{z)dz. 



where F is a closed curve encircling [— c, c] once in counterclockwise directing. Since ^1 is 
analytic at infinity, this integral can be evaluated by means of the residue theorem. Computing 
the residue of ^1 at infinity using Lemma 3.6 then leads to (3.19). □ 

Lemma 3.9. In the multi- critical case a 
(3.21) 



1, r = 1 the following inequalities hold 



Im {Ci,+ i^)-Cl 




) >o. 


X G ( 


-c*,0)U(0,c* 


Im 




) <o. 


X G ( 


-c*,0)U(0,c* 


lm{Cs,Az)-C, 




) <o. 


z £ii 


nw, 


Im(e^_+(z)-e: 




) >o. 


2; G ill 


nw, 


Im {Q^4x)-C 


;,+(^) 


) >o. 


X G 1 


AW, 


Im {^l_ix)-C 


IM 


) <o. 


X G 1 


AW- 
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Proof. In the multi-critical case the (^-function has the simple form 



(3.22) 



w 



i = 1,2,3,4. 



Exploiting this relation we obtain the following algebraic equation for the ^-functions in the 
multi-critical point 



(3.23) 



+ z^ 



0. 



We look for values of z for which this equation has coalescing solutions. Computing the 
discriminant 

z^(256- 27z2) 

we find out that this can only be the case in the points z = and z = itc*. For the case 
z = we clearly have ^*(0) = 0, j = 1,2,3,4. For the case z = itc* only two solutions of 



(3.23) coincide. It follows from Lemma 3.4 that we must have ^J(ibc*) = Q{ic 



We will now prove (3.21). Observe from (3.14)^(3.15) 



= Ci,^ix) = C2,±ix), X G (-c*0) U (0,c*), 

and hence 

(3.24) Im(ei%(x)-C2V(x)) =2ImCiV(x)^0, x € (-c*0) U (0, c* 



where we used that '^*_|_(a;) cannot coincide with ■^2+- Then (3.21) follows from (3.24) and 
Lemma 3.5 by continuity. The other inequalities can be proven likewise. □ 

3.4 Modified A-functions 

The modified A-functions are defined as the following antiderivatives of the modified ^- 
functions 

(3.25) X,{z)= [\j{s)ds, 

Jo 

for j G {1, 4} and z £ I U II, and for j G {2, 3} and z G III U IV. In the other quadrants, 
i.e. for j e {1, 4} and z G III U IV, and for j G {2, 3} and z e lUlI, Xj is defined by 



(3.26) 



X,iz) 



/ S,j{s)ds + i-K. 
Jo 



In both cases, the complete path of integration must stay in the same quadrant. 
The following lemma describes the jump behavior of the A-functions. 



Lemma 3.10. 


The functions Xj, j = 1, 2, 3, 4, 


are analytic within each quadrant 


these jumps 






(3.27) 


Ai^+(x) = Xi-{x) — 2Tri 


X < —c 


(3.28) 


X2,+ {x) = X2-{x) + 2-Ki, 


X < —c 


(3.29) 


^i,+ {x) = Aj_(x), 


x>c, j = l,2. 


(3.30) 


Ai,±(a;) = A2,=F(a;), 


— C < X < c. 


(3.31) 


Ai,+ (2) = Aj-(^), 


z£iR\{0}, i = l,4. 


(3.32) 


X2,±{z) = X^^^{z), 


z G iM\{0}, 


(3.33) 


X'i,±{x) = A4,=f(x), 


X G M \ {0}. 
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Moreover, we have the following symmetry relations 



(3.34) 
(3.35) 
(3.36) 
(3.37) 
(3.38) 
(3.39) 



^2,+ 
A3,+ 

^4,+ 



(x) 


= Ai,- 


(x) 


(x) 


= A2,- 


(x) 


(z) 


= A2,- 


(^) 


(z) 


= A3,- 




(x) 


= A3,- 


(x) 


(x) 


= A4,- 


(x) 



+ iri 

— TTl 



+ TTl 
— TTl 



x G (-c,0) U (0,c), 
x £ (-c,0) U (0,c), 
z G iM\ {0}, 
z G iM\ {0}, 
x G M\{0}, 
X G M\{0}. 



Proof. The analyticity within each quadrant and the equahties ( 3.30 )-( 3.39) are immediate 
from (|3.25|)-(|3.26|) and Lemma Isll Equations (|3.27|)-(|3-29|) are a corollary of (|3.25p-(|3.26 ) 



and Lemma ^ 



□ 



The behavior of the A-functions as z tends to oo follows by combining Lemma 3.6 with 



(3.25)-(3.26). The asymptotics are expressed in terms of functions 9j{z), j = 1,2,3, intro- 
duced in (3.2 ). 



Lemma 3.11. As z ^ oo, we have that 

,2 



Ai(z) 

A2(^) 



z 



log z + £i + 0{z- 



'\logz + i2 + Dz-^/^ + 0(z-^/3) 



eiiz)+{ 



^ log z + £3 + Dljz-^/^ + 0(z-4/3) 

J log Z + £4 + Z)6<j2z-2/3 + 0(z-4/3) 
' i log Z + £3 + Z)wZ-2/3 ^ 0(^-4/3) 

^ log z + £2 + Dz-^/^ + 0(z-4/3) 

J log Z + £4 + Z)w2z-2/3 + 0(z-4/3) 

I log z + £4 + Da;2z-2/3 + 0(z-^/3) 
i log z + 4 + Duz-^/^ + ©(z-^/^) 

where D G M and j = 1,2,3,4, are constants satisfying 

2 




in lU IV, 
in II, 
in III, 

in I, 

in II U ///, 
in IV, 

in lU II, 
in III U IV, 



^2-^4 



-TTl. 



Note that this behavior resembles very much the behavior of the original A-functions in 
Lem. 4.14]. This will be essential in Section 4.3 to perform the second transformation in 
the steepest descent analysis. In fact, the modified ^-functions were precisely defined in order 
to achieve this. 

We will also need the behavior of the A-functions around zero. 
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Lemma 3.12. In a neighborhood of zero, we have 



(3.40) \i{z) = z'^K{z) + 



X2{z) = z'^K{z) + < 



\^{z) = z^K{z) + <^ 



M{z) = z^K{z) + 



F{z)z^/'^ + G{z)z + H{z)z^l^ 

F{z){-zYl^ - G{Z)Z + H{z){-zfl^ + TTi 
-F(Z)(-Z)V2 _ G{Z)Z - H{z){-zf/^ + TTi 

F(z)(-z)i/2 _ G{z)z + + vri 

-F(z)zV2 + G{z)z - H{z)z^l'^ 
F{z)z^/^ + G{z)z + H{z)z^/^ 
F(z)(-z)i/2 _ G{z)z + F(z)(-z)3/2 + vri 

-F(z)(-z)l/2 - G{Z)Z - H{z){-zf/^ + TTi 

F{z)z^l^ + G{z)z + H{z)z^l^ 
-F{z)z^/^ + G{z)z - H{z)zy^ 

-F{z)z^/^ + G{z)z - H{z)z^/^ 
-F(z)(-z)i/2 - G{z)z - H{z){-zfl'^ + vri 



in / U /I, 
in /// U /y, 

in /, 
in /I, 
in ///, 
in /y, 

in /, 
in /I, 
in ///, 
in IV, 

in lU II, 
in III U IV. 



where 


F, G, H, K 


are even analyt 


(3.41) 


m 




(3.42) 


G(0) 




(3.43) 


H{0) 


= e-/4^-5/4 fl. 



_!_ 

47 



127 



F*(0) 
G*(0) 
F*(0) 



0, 
0, 

2„7ri/4 



Proof. The function Ai(2:^) is clearly analytic in the first quadrant. In fact, following the 
analytic continuation in counterclockwise direction we see that Ai(z^) can be extended to a 
function that is analytic in a neighborhood of the origin and vanishes at the origin. This 
means that there exists even analytic functions F, G, H, K such that 



(3.44) 



Ai(z) = z^/'^F{z) + zG{z) + z^''^H{z) + z'^K{2 



for z in / U //. Hence we have the top equality in ( |3.40[ ). By analytic continuation we also 
obtain the expressions for \2{z) in ///, \z{z) in IV , \ii{z) in I U //, A3 in ///, and finally 
A2 in IV . Hence we found half of the expressions in the lemma. 

To get the other half of the expressions, we note that we must have \j{z) = Xj{—z) + cq 
for some constant cq, depending on j S {1, 2, 3, 4} and the quadrant to which z belongs. This 
follows from (3.16) and the fact that \j is an anti-derivative of ^j. It is not difficult to evaluate 
the constants by taking the limit z — )• 0. This gives the other half of the expressions. 

Finally, the values F(0), G(0), and -ff(O) can be found by comparing the asymptotics for 
^j{z) as 2; — )• as given in Lemma 



3.5 



□ 



As n — )• 00, the functions Aj converge to A*, which are the A-functions associated with the 
critical parameters a = —1 and r = 1. In the following lemma we give an estimate on the 
rate of convergence, which will be useful later on. 
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Lemma 3.13. There exist constants ci,C2 > such that 
\Xi{z)-XUz)\ < cin-^/^\z\^/\ 

\Xj{z) - X*{z)\ <C2n-^/^max{\z\^/^,\z\^/^), j = 2,3,4:, 
hold for z £ C and n sufficiently large. 

Proof. We start by giving an estimate for the functions ^j. If we stay away from the singu- 



larities z = 0,c« it is clear from (3.7), (3.9), and (3.12) that 



|^,-(.)-e*(z)|=0(n-i/3) 



uniformly as n — )• oo. Combining this with the asymptotic behavior of (,j{z) as z — )• and 
z ^ oo given in Lemma 3.5 and Lemma[3.6[ we find there exists constant ci and C2 such that 



Mz)-ei{z)\<c,n-yM''^', 



and 



\^j{z) - a{z)\ < C2n-^/^ max{\z 



□ 



for z € C \ {0}. The statement now follows integrating these expressions. 
3.5 Modified equilibrium problem 

Before we continue, we end this section with a discussion on the interpretation of the modified 
functions. This discussion is only to clarify the motivation behind the definition and is not 
necessary for the upcoming analysis. 

Inspired by ^24, eq. (4.27)], we define measures by 



(3.45) 





1 




^(-) 


1 






_ 1 


27ri 



2^ (^'l, 



X G [-c,c], 
z G zM, 
X € M. 



Although (3.45) is exactly the same as [24] eq. (4.27)], the measures here are different because 
our definition of ^ is different. Only in the critical case a = —1 and r = 1, the definitions 
agree. In that situation, the vector of measures {fJ-i, ^2, fJ-s) is the unique solution to the 
equilibrium problem in Section 2.2 See Figure 6] for a plot of these critical measures. For 
other values of a and r, the measures defined by ( 3.45[ ) do not solve that equilibrium problem, 
but a modification of it that we will briefiy discuss. 

Let us first make some observations. It can be checked that the total masses of 
j = 1,2,3, are given by 



/ui([-c,c]) = 1, 



H2m 



and H3{R) 



1 

3' 



e.g. for /ii this follows from Lemma 3.8 So the measures still satisfy conditions (b), (c), and 
(d) of the original equilibrium problem. On the other hand, these measures are not necessarily 
positive. In a neighborhood of zero, and only there, they can have negative densities. The 
measure 112 can also break the constraint locally around zero, violating condition (e) of the 



24 




Figure 6: Densities of the equilibrium measures ^2 (solid), and (dashed) for critical 
values a = —1, r = 1. 



original equilibrium problem. Moreover, generically the densities in (3.45) blow up around 
zero, to either plus or minus infinity. 

Without proof we claim that (//i, //2, Aia) solves the following modified equilibrium problem. 



We seek to minimize the energy functional (2.5) under the conditions (a), (b), (c), and (d). 



but in a small neighborhood of the origin we allow the measures to be negative and ^2 to 
break the constraint a2- These modifications lead to violations of the original equilibrium 
problem, and as a result they will cause trouble when opening the lenses in the steepest 
descent analysis. However, these problems only take place near the origin, around which we 
are going to make a special parametrix anyway. 



4 Riemann-Hilbert steepest descent analysis 



In this section we prove Theorem |2.3[ The proof starts with a characterization of the kernel 
Kn in terms of a 4 x 4 RH problem. We then perform a Deift/Zhou steepest descent analysis 
on the RH problem to obtain the asymptotic behavior near the critical point. 

Many definitions and transformations are analogous to the ones in [24j. For this reason 
we will allow us to be brief at some points and simply refer to for more details. The 
most important differences between the analysis here and the one given in [23], are the use of 
the modified A-functions from the previous section and the special parametrix near the origin 



based on RH problem 2.1 



4.1 Riemann-Hilbert problem for Y 



We start by expressing the kernel Kn in (2.2) in terms of the solution to a RH problem. This 



RH problem first appeared in [38] to characterize the biorthogonal polynomials in (2.1). 



Lemma 4.1. Suppose n = mod 3. There exists a unique 4x4 matrix-valued function 
Y{z), also depending on the parameters a and t, satisfying 

(1) Y{z) is analytic for 2; G C \ M; 
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(2) for X G M the function Y makes the jump 



YUx) = Y-(x) 



(\ WQ^n{x) Wx^n[x) W2,n{x)\ 

1 

1 

Vo 1 y 



where 



-nV{x) I yjg-n{W{y)-Txy) 



dy, j = 0,1,2; 



(3) as z ^ oo, we have 

Y{z) = {I + 0{l/z)) diag (^z"", z'""/^, z'""/^, z-"/^) 



Moreover, the kernel Kn is given by 



(4.1) 



Kn{x,y) 



1 



2'Ki{x — y) 



(0 WQ^n{y) Wi^niy) W2,niy))Y^^{y)Y+{x) 



n\ 






Proof. For details see [211 Sec. 1.7] and the references therein. 



□ 



In the rest of this section we perform a Deift/Zhou steepest descent analysis on Y, i.e. a 
series of explicit transformations 

Y^X^U^T^S^R, 

such that the final RH problem for R is simple in the sense that all jump matrices tend to 
the identity matrix as n — t- cxd, both uniformly and in L^-sense, and the behavior at infinity 
is normalized. Then, by standard considerations, we obtain an asymptotic expansion for R 
as n — )• oo. By tracing back the transformations, we then obtain the asymptotic behavior of 
the kernel Kn- 

The successive transformations are: a preliminary transformation (exactly the same as in 
|24j ) . normalization at infinity by means of the modified A- functions that were established 
in Section [3| opening of unbounded lenses, opening of local lenses, and finally the matching 
with local and global parametrices. 

4.2 First transformation: Y X 

The first transformation Y X is identical to the first transformation in [24, Sec. 5]. For 
that reason, we will not go into details here but only state the new RH problem for X and 
the formula that expresses the kernel in terms of the solution to this RH problem. 



Lemma 4.2. Suppose that n = mod 3. Let 6j be as defined in Section 3.1. There exists 



a unique 4x4 matrix-valued function X{z), also depending on the parameters a < and r, 
satisfying 



(1) X{z) is analytic for z eC\{RUil 



26 



(2) X^{z) = X^{z)Jx, for z G M U iM, where the matrices Jx o,re given as follows 
• On the real line we have for x £ {—oo, —x*{a)] U [x*{a),oo) 

\ 

n 1 n n 

Jx 



l\ g-n(i'V2-ei(x)) Q Q 

1 

e"*^^2.+-f3.+)(a;) X 

VO gn(e3.+-e2,+)(x)y 



and for X £ {—x* (a), x* (a)) 

/I g-n{x'2/2-0i(x)) Q 

Jx = 








1 






\ 



1 g-n(92-e3)(x) 

1 / 



5*66 (2.6) for the definition of x*{a). 
On the imaginary axis we have 



Jx 



(\ 
1 

VO 













jn(ei,_-e2,-)(z) Q 
1/ 



z G 



(3) z — )• oo in the j-th quadrant, we have 

X{z) = (/ + 0(z-2/3)) diag(z",z-'^,z-t,z- 



where the matrices Aj are given by 
(4.2) 



A, 



(4.3) 



A. 



^/3 



%/3 



-1 UJ 
-1 1 1 
-1 a;^ a; 



-1 LO 

1 -1 1 

-1 ^2 



^2 



^4 



V3 



V3 





















LO 


1 






1 


■) 




1 




1 




LO 




-1 


1 




— W 





and w = e2'^*/3. 
(4) X(z) remains bounded as z ^ 0. 



Moreover, the kernel (4.1) is given by 



(4,4) Knix.y) = 2,i(l_y-f (" c-"<»V2-ei(i/)) o) X;'(!/)X+(i) 








Proof The uniqueness for X is standard and the solution can be constructed out of Y by the 
transformation Y X in ^ Sec. 5]. For (lO) see [24, eq. (9.28)]. □ 
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4.3 Second transformation: X U 



In this transformation the modified A-functions of Section 3 come into play to normahze the 
behavior at infinity. As in [241 sq. (6.1)], but now with the modified A-functions, we define 



(4.5) U{z) = {I- nDE2,4)e-'^^^'^X{z)e'^^^'^ 
where 



-nz^/2 







new . ^GC\(MUiM), 



L{z) = { 



and 



diag(^i,^2,^3,^4) for 2 in I, 
diag(^i,^3,^2,^4) for z in //, 
diag(£i,^4,£2,^3) for z in I/I, ' 
diag(£i,^2,^4,^3) for z in IF, 
A{z) = diag {Xi{z), X2iz), Xs{z), X^z)) , 

Q{z) = dmg{ei{z),e2{z),es{z)). 



Here, D is the constant from Lemma 3.11 and S2,4 is the 4x4 matrix with a 1 in the (2,4) 
entry and zeros in all other entries. 

Then U satisfies the following RH problem. 



Lemma 4.3. Suppose that n = mod 3. Then U as defined in (4.5) has the following 
properties 

(1) U{z) is analytic for z £C\{RU iR); 

(2) C/+(z) = U-{z)Jij, for z € MU iM, where the matrices Jjj are specified as follows. 

• On the real line we have for x G (— c, c) 



Ju{x) 



\ 



jn(Ai,+-Ai,_){a;) 







1 

gn(A2,+-A2,-)(x) Q Q 

e'^(^4,+-A4,-)(a;) J 



and for x G (— oo, — c] U [c, oo) 
Ju{x) 



/\ e"('^2,+-^i,-)(a;) Q \ 

1 

^Ti{\z^+-\z-){x) Y 

yO e»i(A4,+ (a;)-A4,_)(x) j 



On the imaginary axis z £ iM. we have 







/I 

1 g™(A3,+ -A3,_)(^) Q 

Vo 1/ 
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iR 



/I 




Vo 




1 

l(A2,+ -A2,_) 












g"(-'*3,+ --''3,-) 2. 

e"(->'4,+-A4,-) 



/l g"(>'2. + -Al,_) 

' 1 



vo 






o»(->^3.+ -A3,_) 








1 

„n(A4.+ -A4 _) 



Figure 7: Jump contour for U, consisting of the axes M and iR, shown together with their 
respective jumps. 



(3) As z ^ oo in the j-th quadrant, we have 

U{z) = (l + 0{z-^/^)) diag(l, z^/\ 1, 



1 
A, 



with Aj, j = 1,2,3,4, as in (4.2)-(4.3) 



The jump contour and the jump matrices are shown in Figure 

Proof. The proof follows by straightforward computations and is similar to the analysis in 
Sec. 6]. □ 



Note that by Lemma 3.10 we have that the diagonal entries of the jump matrices Ju are 
either 1 or highly oscillating. In the next transformations we will open lenses to replace the 
jump matrices containing oscillating entries by jumps for which the jump matrices are either 
constant or exponentially close to the identity as n — )• oo. 



4.4 Third transformation U ^ T: opening of unbounded lenses 

In this section we deal with the oscillatory entries of the jump matrix on iM and in the lower 
right block of the jump matrix on M. We postpone the handling of the jump in the upper left 
block of the jump matrix on the compact set [— c, c] until the next section. 

We open unbounded lenses around iM and M that we denote with L2 and L3 respectively. 
The lips of both lenses emanate from zero and do not intersect elsewhere. The lenses are 
symmetric with respect to both the real and imaginary axes. The unbounded lenses are 
shown in Figure [8| We will be more precise on their definition later on. The transformation 
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Figure 8: The jump contour after opening unbounded lenses around the real and imaginary 
axes. 



?7 I— 7- T and the RH problem for T are essentially the same as in |24^ Sec. 7.1], with the only 
difference that here we use the modified A- functions. 

The definition of T is based on the diagonalization of blocks in Jjj, as explained in 
Sec. 7.1], and given by 



f/l 





(4.6) 



U X < 




1 









1 

1 



[\0 





_g"{A3-A2) 

1 




gn(A3-A2) 

1 








,v 








in the left part of the lens 
around iR, 



in the right part of the lens 
around iM, 



(4.7) 



U X < 



1 
1 










1 




\0 







1 

gn(A3-A4) 








-e'"^"^ 1 J 

o\ 



1 

gn(A3-A4) 



V 



in the upper part of the lens 
around M, 



in the lower part of the lens 
around M, 
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and 

(4.8) T = U elsewhere. 

Lemma 4.4. Suppose that n = mod 3. Then T has the following properties 

(1) T is analytic in C\ Tit, where St is the contour consisting of the real and imaginary 
axes and the lips of the unbounded lenses ; as shown in Figure ^ 



(2) T+ = T^Jt on J^T where Jt is specified in (4. 9) -(4. 12); 



(3) T{z) = (/ + 0(z-V3)) ^-^g A ^ ^1/3^ ^-1 



1/3 



1 

A, 



as z —7- oo in the j-th quadrant. 



The jump matrices have the following structure 



(4.9) 



Jt = < 




on M and on the lips of the lens 
around M, 

on iR and on the lips of the lens 
around iM, 



where the 2x2 blocks are given by 



(4.10) 



grt(Ai,+-Ai,_) gn(A2,+-Ai,_) 
g«(A2,+-A2,-) 



on 



on the lips of the lens 
around M, 



(4.11) (Jt)2 = < 



I g"(A3-A2) 

1 



on lii 



on the lips of the lens around iM, 



(4.12) (Jt)3 = < 



-1 o; 

1 

gn(A3-A4) ^ 



on 



on the lips of the lens around 



Proof. The proof follows by straightforward computations and is similar to |24l Sec. 7.1- 
7.3]. □ 



The precise location of the lips of the lens is given in the following lemma. 
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Lemma 4.5. We can and do choose the unbounded lenses L2 around iR and L3 around M 
such that 



(4.13) 
(4.14) 



Re(A3 — A2) < — C3 min 
Re(A3 — A4) < — C3 min 



^|3/2 I |4/3 



for z on the lips of L2 , 
for z on the lips of L3 , 



for some constant C3 > 0. 



Proof We only prove the statement for the unbounded lens L3 around the real line. The 
other estimate follows similarly. 



For X e M\{0} we know, see (|3^-(|3^ and (|3^-(|3^, that Re(A^^(x) - A^_i(x)) 



0. Then it is clear from the Cauchy-Riemann equations and Lemma |3.9| that the lens can be 
opened such that 

Re(A^ - A^) < 0, zeL^X M. 



The fact that also the stronger estimate (4.14) holds is due to the asymptotics in Lemma 
3ITI Lemma [3I2I and (IsIl-dS^. □ 



Concluding, we have chosen the lips of the lenses such that, for the critical case a = —1 
and T = 1, the off-diagonal entries for the jump matrices Jt on the lips L2 and L3 are 
exponentially small as n tends to infinity. However, we also like this to be true in the non- 
critical situation. This is too much to ask for, but it does hold away from the origin as we 
will show. 

We exclude a shrinking disk 



D (0,n-^ 

centered at with radius > and 
(4.15) 



z £ C \ \z\ < n 



-S 



6 3 



For the purpose of opening lenses, it will be crucial that the excluded disk D (O, n~^) is large 
enough and for this reason we need S < ^. However, when we construct the local parametrix 
m D (0,n-^) we will need that the disk is shrinking with a sufficiently large rate and therefore 
we also need S > ^. The precise value g of the lower bound will be motivated later on, however 
the inequality S < ^ yields the following estimate. 



Lemma 4.6. There exists a constant C4 > such that 



(4.16) 



(4.17) 



Re(A3(^) - X2{z)) < -C4n-i/2j^ax(l, |z|^/3). 



Re(A3(2:) — Xiiz)) < —c^n 



-1/2 



maxfl, \z 



4/3 ^ 



for z on the lips 
of L2 but outside 
D (0,n-^ 

for z on the lips 
of L3 but outside 
D (o,n-^ 



for sufficiently large n. 
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Proof. We prove (4.17). Choose z on the Hps of L3 but outside D (0,n '^). First assume 
\z\ > 1. Combining Lemma 3.13 with Lemma 4.5 gives 



Hence, there exists a constant ci > such that for sufficiently large n 

Re(A3(z) - X^iz)) < -n-^/^ci\z\^/^ 
which proves (4.17) for \z\ > 1. Next suppose that n"^ < \z\ < 1. Then we obtain 

ReiXsiz) - Xi{z)) < (2c2n-^^^ - C3|z|) . 

Now using that \z\ > and 5 < 1/3 we find that 

Re(A3(z) - Xi{z)) < -C2n~^\z\^/^, 

for a certain constant C2 > and large n. Applying the inequality \z\ > once again 
finishes the proof of (4.17). The proof of (4.16) goes analogously. □ 



4.5 Fourth transformation T ^ S: opening of local lenses 

In this section we open the local lens Li around [— c, c]. We let four lips emanate from zero, 
coincide with the lips of the unbounded lens L3 in a small fixed disk D(0,po) around the 
origin, and end in ibc going through the upper/lower half plane. This is shown in Figure |9j 
We will be more precise on the definition of Li later. 
We define 



(4.18) 



T X I 



and 



if 1 o\ 

_gn(Ai-A2) 10 

10 

1 

1 o\ 

en(Ai-A2) 100 

10 

ly 

S = T, 



in the upper part of the lenses around 

[-c,0] and [0,c], 



in the lower part of the lenses around 
[— c, 0] and [0, c], 



elsewhere. 



Then S is defined and analytic in C \ S5 where S5 is the contour consisting of the previous 
contour and the lips of the lenses around [— c, 0] and [0,c]. 
Then S satisfies the following RH problem. 

Lemma 4.7. Suppose that n = mod 3. Then S has the following properties. 

(1) S is analytic in C \ S5; 

(2) 5+ = S-Js on Tis, where Js will be specified in (4.19)"(4.21 ); 
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Figure 9: Jump contour S5 for S consisting of the lips of the local lenses around [— c, c], the 
lips of the global lenses and the axes. Note that the lips of the local lens coincide with the 
lips L3 of the unbounded lens in a small fixed disk D{0,po) around zero (dashed). 



(3) S{z) = (/ + diag (l, z^^^, 1, z"^/^) ^ J as z ^ 00 in the jth quadrant. 

The jump matrices Js are as follows. 

• For z E iR. and z on the lips of L2 we have Js = Jt- 

• For z G M and z on the lips of Li and L3 the jump matrices have the structure 



(4.19) 



Js = 



where the 2x2 blocks are given by 



{Jsh 

{Js)3j ' 



(4.20) {Js)i = { 



-1 0^ 

1 

gn(Ai-A2) ^ 

{Jt)i, 



z e [-c,c\, 



on the lips of Li, 



on M\ [— c, c] and the remaining part of the lips 
ofLa. 
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The block {Js)3 is given by 



(4.21) 



(Jsh 



I2 on the parts of the lips of Li that do not coin- 
cide with the lips of L^, 
{Jt)3 on M and the lips of L^. 



Proof Again the proof follows by straightforward computations and is similar to [24, Sec. 
7.4-7.5]. □ 

We choose the lips of the lens Li such that we have pointwise exponential decay for the 
jumps as n — )• cxD in the critical situation a = —1 and r = 1. 

Lemma 4.8. Let e > be a small constant. We can and do choose the lenses Li around 
[— c*,c*] such that their lips coincide with the lips of L3 inside a disk D{0,po), and such that 

(4.22) Re{Xl{z) - Xliz)) < -c^lzf^^, 

for z S Li \ (R U D{—c*, e) U D{c* , e)) and some constant C5 > 0. 



Proof. The proof is similar to that of Lemma 4.5 



For X e (-C*, 0) U (0, c*) we know, see ([3^-([3^ and ([3^-([3^, that Re(A^_i(x) - 
A2*,±(x)) = 0. Then it is clear from the Cauchy-Riemann equations and Lemma 3.9 that the 
lens can be opened such that 

Re(A^(z) - A^(z)) < 0, zeLi\R. 



Using Lemma 3.12 we see that we can choose Li such that we have (4.22). The fact that 



we can let the lips of the local lens coincide with the lips of L3 in a neighborhood D(0,po) 



of zero is also due to the asymptotics in Lemma 3.12, Here we have to choose po sufficiently 
small. □ 

The following lemma ensures that the jumps on the lenses decay uniformly outside the 
disks D{c*,e), D{—c*,e), and D(^0,n~^) around c* , — c*, and 0, even for the non-critical 
situation. 

Lemma 4.9. There is a constant cq > such that for every sufficiently large n we have 
(4.23) Re(Ai(z) - A2(z)) < -cen^^/^ 

for z on the parts of the lips of Li lying outside the disks D(c* , e), D{—c*, e), and D (O, n"*^) . 



Proof. The proof is similar to that of Lemma 4.6 Because of Lemma 3.13| and Lemma |4.8| 
there exists a constant ci > such that for n large enough 



Re(Ai(z) - X2{z)) < ( 



Now we use that |2;| > n ''to get 

Re(Ai(z) - X2{z)) < -d2\z\^/^n-^ < -d2n-^/^, 
for a certain constant C2 > and sufficiently large n. This proves the lemma. 



□ 
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Now let us stop for a moment and discuss the asymptotic behavior of the jump matrices 
Js as n ^ oo. From (14. 201), (14. 211), (|4.10|)-(|4.12|), and Lemmas |46| and [491 it is clear that the 



jump matrices on the lips of the lenses Li, L2 and -L3 are exponentially close to the identity. 
The same holds true for the jump matrices on (—00, — c) and (c, 00). We leave it to the reader 
to check this. In fact, if we stay away from the points 0, itc these approximations are uniform 
and in sense. This implies that these jumps can be ignored in the asymptotic analysis. The 
final steps in the steepest descent therefore are the following. First, we ignore all exponentially 
small entries in the jump matrices and solve the RH problem that we thus obtain. This global 
parametrix will be a good approximation to the solution to the RH problem for S away 
from the endpoints 0, ibc. Near these endpoints we define local parametrices separately. By 
combining the local and global parametrices we thus construct a function that is a good 
approximation to S everywhere. 

4.6 Global parametrix p(°°) 

If we ignore all exponentially decaying entries of J5 we obtain the following RH problem: 

RH problem 4.10. We search for a 4 x 4 matrix-valued function P^°°\z), also depending 
on the parameters a and r, satisfying 

(1) P(°°) is analytic in C \ (R U iR); 



(2) P_ 



{00) 



p(p°) on M U iM, where Jqo is specified in (4.24); 



(3) P(°^)(z) = {l + {z-^)) diag (1, z^/'\ 1, 
rant; 

(4) = O {z~^/^) as z^O. 

The jump matrices are specified as follows 

' / 1 o\ 
-10 
1 

-1 oy 



1 

A, 



as z ^ 00 in the j-th quad- 



(4.24) 



Jo 



on [— c, c]. 





1 








o\ 







1 



















1 









-1 


° 












0\ 










-1 










1 


















V 



on M \ [— c, c], 



on lii 



This RH problem can be explicitly solved as we show in the next lemma. As the explicit 
formulas are not relevant for our purposes, we will give a short proof only. 
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Figure 10: Image of the map w : TZ ^ C in the w plane. The indicated hnes (plain and 
dotted) are the images of the cuts in the Riemann surface TZ under this map. The solid lines 
are the contours where the function k : C — t- C : u; i— ?• k{'w) changes sign. 



Lemma 4.11. RH problem 4-10 has the solution 



(4.25) P(°°) 



w 



T 



w 



7 

711^2 



Yw4 + T 



7 7 

X diag (k^wi), k{w2)-, ^{ws), k^w^)) 



where k{w) = ((w^ +7^)(u;^ — 87^)) and the square root is taken such that rj is analytic 
in the complex w plane with a cut along wi^{[—c,c]) U W2-{i^) U u)3_|_(M), shown in Figure 
ITU and such that 

k{w) = w'^{l + 0{l/w)) 



as w 



00. Here, tui, . . . , ^4 are the functions of z defined by (3.9) and (3.10). 



Proof. The jump behavior follows from Lemma 3.1 and the jumps of k. A careful asymptotic 
analysis of (4.25), based on Lemmas 3.2 and 3.3 respectively leads to conditions (3) and (4) 
in the RH problem. □ 

The global parametrix fails to be a good approximation near z = —c*,0,c*. Near 

these points we need to construct local parametrices, which we will do next. 
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Figure 11: Contour together with the jumps on each ray. 



4.7 Local parametrices around ±c* 

We now briefly discuss the construction of local parametrices and p('^*) in disks 

D{—c*,€) and D{c*,e) around ibc* for sufficiently small e > 0. The idea is to construct 
g^gj^ that it makes exactly the same jumps in the disk D{±c* , e) as S, and such that 

(4.26) P(=^^*)(z) = {I + O (n"i)) P(°^)(z), uniformly for \z =F c*| = e, 

as re — )• oo. This construction can be done in exactly the same way as in [24 1, Sec. 9] using 
Airy functions. As this construction is standard and the explicit formulas are irrelevant to 



the proof of Theorem 2.3 we will omit it here. 



4.8 Model RH problem for the local parametrix around 



In Section 4.9 we will construct the local parametrix around based on the solution M to 



RH problem 2.1 We formulated RH problem 2.1 in the natural form as it appeared also in 
|20j . For our purposes, however, we need to transform M into M as in (2.13). 
Then, -M(C) is the unique solution to the following RH problem. 

Lemma 4.12. The function M(C) has the following properties 

(1) AI is analytic in C\ where is shown in Figure 

(2) M_|_ = M-J^ on Sjg, where is specified in Figure 
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(3) as C ^ oo we have 

M(C) = B-\C)A~^C± (l + OiC^/^)) DiC), ±ImC > 0, 



where i?(C) = B{iC,) and A,B as in RH problem 



2.1 



anc 



diag (e 



V2(0-it(: pm{0+it< p-miO+iK p-miO-iK] 



Ce/, 



(4.27) DiC) = { 

with 
(4.28) 



miO-itC p-miO+iK ^viiO+nC p-V2{0-iK'j 



diag (e^iK)+**'^, e''2(C)-*<C^ ^-miO-^K 



diag (e^iK)+**'^, e~''2{C)-**C^ ^miO-nc, ^ g-»?i(C)+«tC 



r?2(C) = ir2e3-/4(_^)3/2 ^ 2.2e-/^(-C)^/2; 



), C^iv, 



The fractional powers are taken in the same way as in RH problem 2.1 

(4) M is bounded near zero. 

Proof. The jump contours and the jumps are straightforward. For the asymptotic condition, 
we first compute tpi{iC,) and ^2(^C) fo^^ C in different quadrants and ipi,tp2 as in (2.9 )-(2.10 ) 



7?i(C), CG/umu/F, 



Hence we compute that M satisfies 

M{C) = {l + OiC^))B-^iOA~''C±D{C), ±ImC>0, 

as C — ^ CO. Note that = A and B^ = B. Moreover, by putting B^^A^^C± in front, the 
order term turns into 0{(^^^'^). This gives the asymptotic condition. □ 

4.9 Local parametrix P'^^^ around 

In this section we construct a local parametrix P^^^ in the shrinking disk D (O, n~^) where 
we recall that 5 satisfies the inequality in (4.15). We want the local parametrix to be a good 
approximation to S in this disk, i.e. we want 

(CI) is analytic in D {0,n-^) \ S5; 

(C2) makes the same jumps as S onD {0,n-^) n S5; 

(C3) P^^^ matches with the global parametrix p(°°) on the boundary of the disk, i.e. 



p(^\z)=^I + 0{n'-'/^))p'^^\z), \z\ 



n 



as n — 00. 
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We search for such a P^^^ in the form 



(4.29) P(°)(z) = K(z)M(n2/3/(z);ri,r2,n2/3si,n2/3s2,ni/3^) 



X diag ( e 



^nXi(z) n\2(z) „n\3(z) nX4,{z) 



where M is the solution to the model RH problem and / is a conformal map on a neighborhood 
of the origin. The function En is an analytic 4x4 matrix valued function. Defined in this 
way, P(°) can already be seen to satisfy condition (CI) and (C2) above. Indeed, (CI) is 
immediate, possibly after slightly deforming some contours. The jumps follow from Lemma 
(4.29), and the jumps for M indicated in Figure [ll} Finally, we have the conformal map 



3.10 



f{z), the analytic prefactor £"„,, and the parameters rj, Sj and t at our disposal to make P^"-* 
also satisfy the matching condition (C3). It will be necessary to let the parameters depend 
on z and n. To retain the analyticity of P^^^ the dependence on z will have to be analytic. 



We will argue later that also for these complex parameters RH problem 2.1 is still solvable 
as long as n is sufficiently large. 

We now come to the precise definitions, starting with the conformal map /. We recall that 



F, G, and K defined in Lemma 3.12 are even functions that are analytic in a neighborhood 
of the origin. We define 



(4.30) 
(4.31) 

(4.32) 

(4.33) 
(4.34) 



ri{z) 

r2{z) 

s{z) 
t{z) 



^e— {H{z) - {F{z) - Fmz-') , 
Je— {H{z) + {F{z)-F{0))z~'), 



Sliz) = S2{z) 

-iG{z). 



1 „iTi 



/4f(o), 



Note that ri^2, s, and t are analytic in a sufficiently small neighborhood of the origin. In 
the construction of the parametrix, we will use that the functions rj, s, and t are analytic in 
D (O, n~^) and have the asymptotic behavior as n — >• oo indicated in the following lemma. 
Note that rj, s and t depend on n, although we have suppressed that dependence in the 
notation. For this lemma the lower bound in ( 4.15| ) is essential. 



Lemma 4.13. Let ^ < 6 < ^. The functions Tj[z), s{z), and t{z) are analytic for z G 
D (O, n^^) and n sufficiently large. Moreover, we have the following limits 



(4.35) 



lim„^oo rj{n^^z) = 
lim„_i.oo n^/^s(n"^z) 
lim^^oo n^^^t{n~^z) 



56), 



for \z\ < 1 and j = 1, 2. 

Proof. It was discussed above that these functions are analytic in D (O, n~^^ for n sufficiently 
large. As F{z) is even, (|4.3lD-(|4.32D and (|3.43D lead to 



rj{z) = |e-"^/^//(0) + 0(^ 



1 +y/^j)+0{z) 



127 
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as z — )• 0. The limiting behavior for rj then follows recalling that 7,t— )-lasn— )-oo. 
For s{z) we follow the same strategy. Using (4.33) and (3.41) we obtain 



„2/3^1/4 



which in combination with (3.7) and (2.7) yields the statement. 
The limiting behavior for t can be found by realizing that 

lim n^/h{n-^z) = -i lim n'/^Gin-^z) = lim n^/^-f-^/^ f _ i _ 5^4/3 



which follows from (4.34), (3.42), and the first inequality in (4.15). Inserting the limiting 
behavior for 7 and r as given in (3.7) and (2.7), finishes the proof. □ 



From Lemma 4.13 it follows that (4.29) is well-defined (postponing the definition of En for 
a moment). Indeed, it follows from standard arguments that if the solution to the RH problem 
in Lemma [4.12 exists, it depends analytically on the parameters rj,s and t. Combining this 
observation with Lemma 4.13 and Theorem 2.2 we see that for the choice of rj, s, and t that 
we made, the solution to RH problem exists and hence (4.29) is well-defined for sufficiently 
large n. Moreover, we defined rj, s, and t in order to have the following equalities, which are 
the key to obtain (C3). 

Lemma 4.14. We have that 

(4.36) rji (n^/^z;ri{z),n^^^s{z)) ±int{z)z 



{{-zf/^H{z) + (-z)1/2f(z) ± zG{z)) , 
{-{-zf/^H{z) - (-z)i/2i?(^) ± zG{z)) 



z£l, II, 
z G ///, IV, 



and 



(4.37) r]2 (n'^^^z; r2{z),n'^^^s{z)^ ± int{z)z 



n [—z 
n (z 



y^H{z) - z^^Fiz) ± zG{z)) 
^/'^H{z) + z^/^F{z)±zG{z)) , 



z€l, II, 

z e ///, IV, 



for z£ D {0,n~^). 



Proof. We use the definitions of r/j, rj{z), s{z), and t{z) as given in (4.28), and (4.31 )-(4.34). 
The first step in the proof is to realize that 



„37rj/4 3/2 



-rie""''^z"'" + 2e"'-/^s{z)z^'" = iz' 



,1/2 



^I^H{z) - iz^l^F{z) 

_ { {-z f/^H{z) + {-zY/^F{z), z€l, II, 
~ \ -{-zf/'^H{z) - {-z)^/^F{z), z G III, IV, 



which follows from the definitions by a straightforward check. Then (4.36) easily follows 
combining this with (4.28). 
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The statements in (4.37) follow in the same way. Indeed, 



-J2 



(^)e3W4(_^)3/2 ^ 2e"*/S(z)(-z)^/2 = i{-zfl^ {H{z) + z-^F{z)) 



-z^/^H{z) - z^l^F{z) 
z^/^H{z) + z'^/'^F{z), 



zel, II, 

z e ///, IV. 



□ 

Finally, we define the analytic prefactor En by 
(4.38) Eniz) = e-'''^^^'^P^°°\z)C^^AB (n^/^z^ , ±lmz>0. 

Lemma 4.15. The prefactor En is analytic in a neighborhood of the origin. 

Proof Trivially, En is analytic within e ach qu adrant. Moreover, it follows from the behavior 
of p(°°) around zero, stated in Lemma 



4.11 



and ( |2.12D that En{z) = O (2"^/^) as z ^ 0. 
Therefore, it is sufficient to prove that En has no jumps on the real and imaginary axes. 
To this end, we note that by some elementary computations one can check that near the 
origin B{z)~^A-^C± satisfies the same jump conditions as Therefore, En has indeed 

no jumps and hence it is analytic. □ 

We have the following result on the local parametrix P^'^K 

Proposition 4.16. // n = mod 2, then the function P^^^ as defined in (4.29) satisfies 



(C1)-(C3). 

Proof. In the discussion right below the conditions (C1)-(C3) we already argued that the 
conditions (CI) and (C2) are automatically satisfied, so it remains to show that we have 
(C3). 



By inserting (4.36) and (4.37) in (4.27) and using Lemma 3.12 and n = mod 2, we 
obtain 



(4.39) 



D{n^l^f{z)) = e'^^'^(^) diag(e-"^i(^\e-"^2W^g-nA3W^g-nA4(^)^ _ 



Note that we need n = mod 2 because (3.40) and ( 4.36 )-( 4.37) differ by a term vrz at some 



places. Recall that p(°°)(z) = 0{z~^/^) and p(°°)(z)-i = Ojz^/^) as z ^ 0. Then we have 
pW(z)p(°^)(z)-i = P(°°)(z) (1 + 0(n-i/3^-i/2)^ p(°°)(z)-i = / + 0{n^-^l^), 



where the order is uniform for Izl 



n 



□ 



4.10 Final transformation: S R 

We define the final transformation as 



(4.40) 



R{z) 



' S{z) (PW) ^ (z), for z G D (0, n-^) , 
S{z) (p(±c*))"^ (z), for z G D{±c*,e), 
^S{z)P^°°^ ^{z), elsewhere. 
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Then R is analytic in C \ S/j where S^j is the contour consisting of the three circles around 
0, — c*, and c*, and the parts of the real and imaginary axes and the lips of the lenses that 
lie outside these circles. The circles are oriented clockwise. 
Then R satisfies the following RH problem. 

Lemma 4.17. Suppose that n = mod 6. Then 

(1) R is analytic in C \ S/j; 

(2) = R-Jr on TiR with jump matrices 



f p(0) ('p(oo)^ 



Jr 



p(±c*) 



for \z\ = n 
for I z =p c* I 



^ Js (P(°^))^^ elsewhere on ^r; 



(3) as z —7- oo we have 



R{z) = I + 0{l/z) 



Now, all jump matrices are close to the identity matrix as n gets large. On the segments 
of the real and imaginary axes belonging to Tir the jump vanishes. On the shrinking circle 
around the matching condition (C3) in the construction of P^*^^ in Section 



4.9 



yields 



jR = I + 0{n' 



5-1/3 



as n 



oo, 



uniformly for \z\ 



n 



From (4.26) and (4.23) we find 



Jr = I + O in , as n — )■ oo. 



for z on the other bounded parts of the contour, i.e. on the circles around — c* and c* and 



on the remainders of the lips of the bounded lens Li. Moreover, in view of (4.16) and (4.17), 



there exists a constant d > such that on the unbounded parts of the contour T^r the jump 
matrices satisfy 



Jr = I + 



-dnV2max(l,|2|4/a) 



as n 



oo, 



uniformly for z on the unbounded arcs of Hr. 

We have now achieved the goal of the steepest descent analysis for the RH problem for 
Y . R{z) tends to the identity matrix as z — ?■ oo and the jump matrices for R tend to the 
identity matrix as n — )• oo, both uniformly on T,r and in L^(S/j). By standard arguments, 
see [15] and in particular [7| for the case of a moving contour, this leads to the conclusion of 
our steepest descent analysis. 



Lemma 4.18. Let R he as in (4.40) and n = mod 6. Then 



(4.41) 



R{z) = I + 



n 



5-1/3 



as n ^ oo, uniformly for z £ C \ 'Sr. 
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4.11 Proof of Theorem [231 



Proof of Theorem 2.3. We start with the formula (4.4) that expresses Kn in terms of X. 
The idea is to write this expression in terms of R instead of X by applying all respective 
transformations 

X^U^T^S^R 

introduced in the steepest descent analysis. Meanwhile, we let n tend to infinity, so that we 
can exploit the conclusion of the steepest descent analysis (4.41). 



Let X, y G M. First, unfolding the transformation X i— )• [/, given in (|4.5|) yields 

-n{yy2-xy2-\2,+ {y)+Xi,+{x)) 



Kn{x,y) 



2TTi{x — y) 



(0 1 Q)U^\y)U+{x) 



/1\ 





Vo/ 



The opening of global lenses ?7 i— t- T in (4.6)-(4.8) does not effect the above expression. 
Indeed, as x and y are real, they fall outside the global lenses around the imaginary axis. 
Also, the opening of the global lens around the real axis does not affect this expression. Thus, 



-n(j/2/2-zV2-A2,+ (j/)+Ai,+ {x)) 

Kn{x,y) = ^^-(—^ (0 1 0)r+-i(y)T_,(x) 



/1\ 






The opening of the local lens T i— >• 5 in (4.18), however, does have impact: 

-n(y2/2-x2/2) 



(4.42) Kn{x,y) 



27ri(x — y) 



(_e"^i,+(3/) e"^2,+(j/) {)) S^^{y)S+{x) 



/g-nAi,+ (x)\ 
,-nA2,+ (x) 



V / 



if we assume —c<x,y<c. Moreover, for x,y D (0,n , we have 

S{z) = Riz)P^'\z), z = x,y, 



by ( |4.40[ ). It then follows that 

p-nfo2/2-x2/2) 

(4.43) Kn{x,y)-- 



-110 0) 



2TTi[x — y) 

M-i (n2/3y;ri(y),r2(y),n2/3s(y),n2/3s(y),ni/3t(y))i^-i(y)ii-i(y)fl(x)i^„(x) 



M_|. ( n^/^x; ri{x),r2{x),v?/^ s{x),r?/^ s{x),n^/^t{x) 



/1\ 

1 





Now we scale x and y with n such that 

u 

j^2/3 



and y 



n 



2/3' 
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and u, f G M. Then for large n, x and y belong the the disk D (0,n so that (4.43) holds. 
We want to take the limit as n — )• oo. Note that under these conditions 

hm e-"(2^V2-xV2) = 1 



and by (4.35) we have 



1, 



n^/^t(z) ^ -a, 



J = 1,2, 
J = 1,2, 



as n — )■ oo and z = x,y. Furthermore, it follows from (4.41) and Cauchy's formula that 



R-Hy)iRiy) - R{x)) 



(4.44) 



(y - x) / R-Hy){R{w) - I) 
{w — x){w — y) 



O 



21X1 
x-y 



On' 



as n — >• oo where the constant is uniform for u, v in compact subsets of M. 

We estimate the factor En given in (4.38) in the same way. First note that we have 
En{n~'^/'^u) = O {v}/^^ and En{n~'^l'^v) = O {p}/^^ as n — >• oo uniformly for u,v in compact 
subsets. Then we obtain 



E-\y){En{y) - En{x)) 



{y-x) r ( E-\y){En{w) - I) 



27ri 



{w — x){w — y) 



dw 



as n — >• oo where the order term is uniform for u, v in compact subsets. This results in 

lim E;,\y)R-Hy)R+{x)En{x)=I, 



uniformly for u, v in compact subsets of M. Inserting the estimates for R and En into (4.43) 
yields 



lim 



1 



Kn 



U V 



1 



(-1 1 0) 



n— s>oo fi,2/3 \^2/3'j^2/3/ 27ri(u — v) 

X [v- 1, 1, \{a^ - 5b), lia" - 56), -a) M+ {u; 1, 1, ^(a^ - 56), {{a^ - 56), -a) 



/1\ 

1 





This proves the statement. 



□ 



5 Solvability for RH problem 2.1 



The aim of this section is to prove Theorem |2.2[ Note that by a simple rescaling it is sufficient 
to consider the case ri = r2 = 1 only. We will assume this throughout this section. 
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5.1 Lax pair 

Consider first tlie system of equations 



(5.1) 



l^MiC,t) = UiC,t)MiC,t), 
§-^M{C,t) = W{C,t)M{C,t). 



where the matrices U and W are exphcitly defined as follows. We need the Hastings-McLeod 
solution q^v) to the Painleve II equation as defined in (1.2) and (1.3). We will also use the 
Hamiltonian 

(5.2) uia) = {q'{a)f-aq''ia)-q\a). 

Note that u'ia) = -g^(cr). Then we define U and W as 



(5.3) 



and 



(5.4) 



Here 



u{C,t) 





( c 


-2b 


-2id\ 




w{c,t) = 


-26 


-C 2id 










-2if C 


-2h 






\2if 


-2h 


-c ) 




/ t-c 




d 


i 


\ 


-d 




c-t 





i 


i{C — s) + ik 


-i(h + h) 


t + c 


d 


\ -i{h + h) - 


-i{C + s)+ik 


-d 


-{t + c)l 



'd = 2-y^q (22/3(2s-t2)), 

C= -2-1/3^ (22/3(2s-t2)) +s2^ 



(5.5) 



b= ^^f^+dc + td, 



i_dd 

2t at 



f = 2dt^ - c 
h= l^^ + dc-td, 
k 



dc^ -d^ - 2ds, 



c^-d^- s. 



Since q has no poles on the real axis [32], U and W are defined for all t € M. 

Note that the Lax pair is an overdetermined system of equations for ^. However, the 
compatibility condition is satisfied. 



Lemma 5.1. With U and W defined as in (5.3) and (5.4) we have 

d 



(5.6) 



^^w = l-u + [u,w]. 



In other words, U and W solve the compatibility condition for (5.1). 
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Proof. This is intrinsic to the construction of U and W and follows by direct verification. 
Plugging (5.3) and (5.4) into (5.6) it remains to check the following equalities 



c' = 2{h - b)d, 

d' = 4b{t - c) - 2ds + 2dk - 2f, 
= Ah{t + c) + 2ds - 2dk + 2/, 
b-h = 2dt, 

k' = 2{h'^ - 6^), 
h' + b' = -Aft + 2{h -b){k-s). 



These are all valid by virtue of (5.5). 



□ 



The main idea of our proof of Theorem 2.2 is the following. By standard arguments of 



isomonodromic deformation theory, we construct a RH problem out of the Lax pair (5.1) for 



which the jump matrices do not depend on the parameters. It is crucial that the Hastings- 
McLeod solution to the Painleve II equation has no real poles, so that we get a solution for this 
RH problem for every s,t £M. We then prove that this solution also solves the RH problem 



2.1[ which finishes the proof of Theorem 2.2 To this end we need to check that both the jump 
conditions and the asymptotic condition are satisfied. The asymptotic condition follows by 



construction and more specifically the top equation in (5.1). For the jump conditions we use 



the fact that they are independent of t which allows us to simply refer to [20j for the special 
case t = 0. 



5.2 Asymptotics 



We will now show that the top equation of ( 5.1 ) has a unique formal solution near the irregular 
singularity oo that has the asymptotic behavior given in RH problem 2.1 It turns out that 
this equation exhibits a Stokes phenomenon. This means that to find a genuine solution of 
the equation with the same asymptotic behavior as the formal solution, we need to restrict 
ourselves to certain sectors in the complex plane. To ensure existence of such solutions, these 
sectors should be sufficiently small. On the other hand, to get uniqueness the sector needs to 
be sufficiently large. 

Motivated by the asymptotic behavior in RH problem 2.1 we define the function A by 
(5.7) A{C) := B{OA diag (e-^(-^)+*^, e-^(^)-*^ , e^^^-^^+^f, e^(^)-*< 

where A and -B(C) are given by (2.11) and ( |2.12[ ) and 

V.(c) = ic^/2^2<i/2. 

Note that A has branch cuts along the positive and the negative real line. 
Now consider the sets 



where 



V'i(C) 



= {C G c I ReiMO - i^jiO) = 0}, 1 < i < j < 4, 
(-C)'/', MC) = -C'^\ MC) = i-Cf^', MO = 
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Then each ii j is a union of rays that connect the origin with oo, see also Figure 12 These 



rays are called the Stokes lines associated with A{C,)- They separate the complex plane 
into sectors of angle 7r/6. The scheme of dominance is indicated by the numbers between 
successive Stokes lines, e.g. the sequence 4, 1,3,2 in the sector < argC < 7r/6 means that 
ReV'4(C) > Re^/^i(C) > Re^/^3(C) > Re?/'2(C) fo^^ C iii that sector. Note that the Stokes lines on 
M are special in the sense that the order of dominance does not change when crossing these 
lines. This is due to the branch cut of A{C) along the real line. 
We now define sectors il^, /c = 0, . . . 5 in the following way 



(5.8) 



TV TV TtV TV 

C e C \ {0} I + A;^ < argC < ^ + A;^ 



Note that these sectors cover C \ {0} and that there is overlap between Qk and ^k+i- See 
also Figure 12 Note also that the sectors Qj, j = 0,2,3,5 contain the positive or negative 
real line on which ^(C) has a branch cut. 

In order to deal with the technicalities coming from the branch cuts in the definition of 
A we introduce some extra definitions. First, we define ^^(C) a-s the analytic continuation 
of A{C) in the upper half plane to C \ [0,—ioo) and A~{C) as the analytic continuation of 
A{C) in the lower half plane to C \ [0,-100). In the same way as we did before for A we can 
construct Stokes lines and a scheme of dominance. For A'^{C) we find 



= {C e C I Re(V'+(C) - V'+(C)) = 0}, 



1 < f < j < 4, 



where i'tiC) is the analytic continuation of ^j(C) in the upper half plane to C\ [0, —ice). The 



scheme of dominance for ^^(C) is shown in Figure 13 We leave the details for ^^(C) to the 
reader. Note that the sectors Qq, Oi, and 0,2 contain precisely one Stokes ray from ifj for all 
pairs i < j- We refer to such sectors as Stokes sectors associated with A'^. Analogously Q3, 
Q4, and fis contain precisely one Stokes ray from i~j for all pairs i < j. 

We then have the following result on the existence and uniqueness of solutions of the top 



equation in (5.1) with the desired asymptotic behavior. 



Lemma 5.2. Fix t G M and let U be as defined in (5.4). Moreover, let Qj^ be one of the 



sectors defined in (5.8). Then the equation 

dN 



(5.9) 



dC 



(C) = c/(C,t)iV(C), 



has a unique fundamental solution N in the sector 0^ such that N{() has the following 
uniform asymptotics as ^ ^ 00 within this sector 



(5.10) iV(C) 




Here, Ni is given by 



(5.11) 



iVi 





U+(c) 


A; = 0,1,2 




)^-(C) 


A; = 3,4,5 


/ * b 


* id\ 




—b * 


id * 




* if 


* h 




\if * 


-h *) 





where * denotes (for our purposes) unimportant entries. 
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Figure 12: Stokes diagram for the irregular singular point oo of the differential equation (5.9) 



associated with A{C)- The stokes lines ii^ separate the complex plane into sectors of angle 
tt/6. The scheme of dominance is indicated by the numbers between successive Stokes lines, 
e.g. the sequence 4, 1,3,2 in the sector < argC < n/G means that Re'04(C) > > 
Reip3{C) > Re ■02(0 ^oi ( in that sector. 
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Figure 13: Stokes diagram for the irregular singular point oo of the differential equation (5.9) 



associated with ^+(<^). The stokes lines ifj separate the complex plane into sectors of angle 
7r/6. The scheme of dominance is indicated by the numbers between successive Stokes lines. 
The zigzag line indicates the branch cut of A'^{C,). 
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Proof. We prove the lemma for = 0, 1, 2 and leave the other case to the reader. 



Infinity is an irregular singular point of the ODE (5.9). Asymptotics of solutions to equa- 
tions of this type can always be computed following an algorithm presented in the monograph 
|48j . When we consider U as a linear polynomial in (, its leading coefficient matrix is nilpo- 
tent. First we will transform this system in order to obtain a leading coefficient matrix which 
has four different eigenvalues. As soon as we are in that situation we can invoke a general 
result, see e.g. 03| or Th. 1.4], to obtain asymptotics of the transformed system. 

Thus, we start with the transformation 



(5.12) 



where B~^{Q is the analytic continuation of B{Q, see (2.12), in the upper half plane to 
C \ [0, — ioo). We supplement the transformation by the change of variables = z^, where we 
choose z such that — 7r/4 < argz < 37r/4, and call Z{z) = Z{C)- In this way we obtain the 
equivalent equation 



(5.13) 
where 



z 



1 f)7 

^^ = Eiz)Ziz), 



E{z) 



Eq H — El H — ^£^2 H — ^E^. 

Z Z'^ 



The matrices Ej, j = 0,1,2,3 can be explicitly found by a lengthy computation and are 
constant 

Eo = diag {-2i, -2, 2i, 2) = diag (Ai, As, A3, A4) , 



^1 



/ 2t 

-V2id 
-2ic V2id 
V V2id -2ic 



V2id 2ic V2id\ 
-2t V2id 2ic 

2t V2id 

V2id -2t / 



Eo 



( i{s-k) 
- exp(7ri/4)(6 + h) 
s — k 

\-iexp(7ri/4)(6 + h) 



exp(7ri/4)(6 + h) 

k — s 
-i exp(7ri/4)(6 + h) 
i{k — s) 



s — k 
i exp(7ri/4)(& + h) 

i{k — s) 
-exp(7ri/4)(6 + h) 



iexp(7rz/4)(6 + h)\ 

i{k — s) 
exp{m/4:){b + h) 
s-k ) 



/o 





-i/2 


0\ 











i/2 


i/2 











Vo 


-i/2 





0/ 



E:, 



Hence £"0 indeed has four different eigenvalues, which was precisely the goal of the transfor- 
mation. 



The ODE (5.13) is of the correct form to apply [31', Prop. 1.1]. This proposition states 



that (5.13) admits a unique formal fundamental solution 
(5.14) Z,(z) 



AopA3 25-+A2 4r+^i2 



z "e 
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where = /, A^, k = 0,1,2,3, are diagonal, and A3 = Eq. Moreover, all coefficients A^, 
A: = 0, 1, 2, and Z^, k = 1,2,3, .. ., can be determined recursively using the explicit form of 
Ej,j = l,...,4. 

If we unfold the transformation (5.12) this means that (5.9) has a unique formal solution 
for G C \ [0, —zoo) of the form 



(5.15) 



/ \ 3/2 



Here the A^^ are constant matrices, that can be recursively found from the matrices Z^. 



(5.16) 



To determine all coefficients, we rewrite (following [3T]) (5.14) in the equivalent form 

/ 00 

Zf{z)= [Y^YkZ- 



Ya = L 



\k=0 



where all matrices Y^, k = 1,2,3,... are off-diagonal and the diagonal matrix A{z) is given 
by the series 



z^ z^ °° 
A3— + A2 — + Aiz + log Ao(2) + X] 



k=\ 



-k 



Now substituting (5.16) into (5.13) and equating coefficients of the same powers of z yields 
a system of equations for the diagonal matrices A^ and the off-diagonal matrices Yfc. This 
system can be recursively solved, which leads to 



(5.17) 



A' 



3-fc 



diagFs.fc, 



3-k)i 



Xj — \i 



where F^-k are defined by the equations 



(5.18) F3-k = Ek+Y,iFk-mYm-Y^A3-k+m)-{k-3)Yk-3, k = 0,1,2,.... 



m=l 



Here we have put Ek = 0, k > 3, and Yk = 0, k < 0. In this way we find e.g. 

A3 = Eo, 

A2 = diag{2t,-2t,2t,-2t) , 
Ai = diag {2is, —2s, 2is, 2s) , 
Ao = 0, 



where we used the last line of (5.5). 

By continuing this approach we can in principle obtain all coefficients in (5.15). We claim 
that all coefficients iV{2'+i), / = 0,1,2,... vanish. To see this, observe that the matrices 
AEkA~^ are block diagonal (with blocks of size 2x2) when k is odd and block anti-diagonal 
for even k. It then follows by ( 5.17| )-(5.18) that AYkA^^ and AAkA~^ are block diagonal for 
even k and block anti-diagonal for odd k. In the context of (5.14) this means that also Zk is 
block diagonal for even k and block anti-diagonal for odd k. Taking this into account when 
unfolding the transformation (5.12) proves the claim. 
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Now we proved that (5.9) admits a formal solutfon of the correct form. The fact that in 
each sector fi^, k = 0,1,2 this formal solution actually has the asymptotic meaning stated in 
the lemma follows from the choice of the sectors O^, see \3T\ Th. 1.4]. 



It remains to prove that iV(2) = Ni as in (|5.1ip. This follows from this approach by 



detailed but straightforward calculations in which we also use (5.5). Alternatively, (5.11) can 
be checked by substituting the asymptotic formula (5.10) into (5.9). □ 



5.3 Proof of Theorem Q 



We are now ready to prove Theorem 2.2 



Proof of Theorem \2.S\ The proof follows by standard isomonodromy deformation arguments 
and the results in [20j for the case t = 0. For more details on the general theory we refer to 

m- 

First we remark that without loss of generality we may assume that the angles in RH 
problem 2.1 coincide fi = f2 = 7r/3 such that T,m = u|^Qe^'^*/^[0, cxd) with constant jump 



matrices on each of these rays. The rays separate six sectors 



A, 



From Lemma 5.2 



{CgC\{0} I fcf <argC< (A; + l)f}. 
we learn that in each sector Qj^ we can find a unique solution ^fc(C) t) to 



the top equation of (|5.1|) that has the asymptotics given in (|5.10|) as C 
that if we let C 



oo within Qi.. Note 



oo in the smaller sector this coincides with the (uniform) asymptotics in 
the RH problem. Moreover these functions can be analytically continued on the full complex 



plane. The solutions associated with different sectors differ by a matrix Sk{t) 
k = 0, . . . ,5, that is independent of Finally, we define a function M{(; t) by 

M(C;t) = ^fc(C;i), CGAfc. 



Summarizing, we have constructed a function M(C; t) with a jump contour on T,m that satisfies 
the asymptotic condition as C — ^ oo. It remains to show that the jump matrices Sk{t) on Sjvf 
coincide with the jump matrices in the RH problem. 

For t = this follows from the results in [20]. We know from [2^ that RH problem 2.1 
has a solution M{() for t = and that this solution satisfies (5.9). Now define ^k{C) &s the 



analytic continuation of M{(^) in the sector A^ to the full complex plane. It can then be 
checked that ^a:(C) uniformly satisfies the asymptotic condition stated in (5.10) as — )• oo 
within flk, and as a result we have ^'a:(C) = ^A:(C; 0) and hence M{() = M(C; 0). This means 
that for t = 0, we indeed have that 5*^(0) are the jump matrices in the RH problem 2.1 



By using standard arguments from the theory of isomonodromy deformations we will now 
show that the Stokes matrices Sk are independent of t and that finishes the proof for general 
f G M. To this end we note that from (5.6) it follows that = d^^/dt — W^k also solves 
the top equation of (5.1). This implies that there exists a matrix C{t), not depending on 
such that ^'fe 



^fcC(t) and hence 
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We then show that C{t) vanishes. We assume here k = 0,1,2. The adaptations to cover 
fc = 3, 4, 5 are straightforward. Computing the asymptotic behavior for (" — )■ oo of the right- 



hand side, using also (5.11), yields 



X diag ( e 



-H-O+K ^-HO'K ^'H-O+K ^'m-K 



as C — ^ oo in fi/^ 



where (again) the superscript + indicates the analytic continuation of the function in the 
upper half plane to C \ [0, —ioo). Letting ^ — t- oo within 0,^ this already shows that C{t) has 



entry of C{t) is zero it is crucial that 
12 Therefore contains a subsector 



zero diagonal. To see that also the (i, j) off-diagonal 
Qk has a Stokes ray in ifj in its interior, see Figure 

in which Ke{ipj~{C) — ''PtiO) < 0- Letting C approach infinity in this sector proves that 
C{t) 



0. Hence C = 0. 



We thus proved that also solves the bottom equation of (5.1). Therefore the Stokes 
matrices are independent of t. Indeed 



(5.19) 



and hence dS^/dt = 0. This finishes the proof. 



□ 



6 A double scaling limit for Kc 



In this section we prove Theorem 2.6, The proof is again based on a Deift/Zhou steepest 



descent analysis on the RH problem 2.1 for M. Since the analysis is much simpler than the 



one in Section 4, we will allow ourselves to be brief. We will give all the transformations, but 
we will be short about the proofs. 

The steepest descent analysis consists of a sequence of transformations 



In the end, the asymptotics for M^^) (as a — )• oo) can be easily found. The proof of the 
statement then follows by updating the kernel with each transformation. The core of the 
proof is that at a certain point, we need to construct a local parametrix, which is done using 



RH problem 2.5 



6.1 First transformation: M t-^ M^^^ 



In Theorem 2.6 the solution M(z) = M{z;s,t) to RH problem 2.1 appears in the critical 
kernel for the values of parameters 



n =r2 = 1, 
t = -a{l-^) 

S = Sl = S2 = 



2 I 1 

2 



a 

2 ' 



where cr G M is fixed. 



54 



We define the first transformation as 
(6.1) M«(z;a) = diag (a^/^, a^/^, a'^/^, a-^/^) M {a^z; (l - ^) 



Then from RH problem 2.1 for M, we straightforwardly check that M*^^) satisfies the following 
RH problem. 

Lemma 6.1. The function M^^) defined in (6.1) has the following properties 

(1) M(i) is analytic for ( £C\ T.m; 

(2) M^^\z) = M^}\z)Jk, forzeTk,k = Q,..., 9; 

(3) z — 7- oo we have 

M^^\z- a) = {l + (z-i)) diag ({-z)-^l\ z-^l\ ^-zfl\ ^^4^ ^ 

xdiag ^ffl\-^{^)^^^-'^I^M ^e^ 

where 

4>{z) = \z^l' + z^/^- 

(4) M(^)(2;) is bounded near z = 0. 

6.2 Second transformation: M^^^ ^ M^^) 

In the second transformation we deform the contours. We shift the rays Fi and Fg of the 
contour Tim to the right by one, while shifting r4 and Fg to the left by one. More precisely, 
we introduce the new rays 

Fi,9 = Fi,9 + 1, and F4,6 = F4,6 - 1. 

Moreover, for A; = 0, 2, 3, 5, 7, 8 we set F^ = F^ and define 

9 _ 

k=0 

We define the function M^^^ as follows. For k = 1, 4, 6, 9 it is defined by 

(6.2) M(^)w=^';;;'--'''\ '^=''''' 

• ' ' ' \mW(2)J,-' fc = 4,9, 

for z in the region enclosed by F^, F^ and M. In the other regions we simply set M^^) = m(^). 
Lemma 6.2. The function M^'^^ as defined (6.2) has the following properties 

(1) M(2) is analytic /or C G C \ ^^(2); 

(2) Mf (z) = Mi')(z)Jjv,(2), for z e 

(3) M(^)(z) /las the same asymptotics as M^^^ as z ^ oo; 

(4) M^'^\z) is bounded near z = 0. 

See Figure [T^] for the jump matrices Jf^^cz) ■ Note that we reversed the orientation of some 
rays. 
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11/ 


\0 





Figure 14: The jump contour and jump matrices for M^'^\z;a). 
6.3 Third transformation: M^^) M^^) 

In the third transformation we use ^-functions to normahze the behavior at infinity. The 
^-functions are defined by 

9i{z) = - zf/-' - (1 - a/a^)z, g^iz) = + zf^ + (1 - a/a^)z, 

(6.3) 

g^{z) = \{l- zfl^ - (1 - ala^)z, g,{z) = \{l + zf/^ + (1 - a/a^)z. 

These functions satisfy the following properties; 

• gi and 53 are analytic on C \ [1, +00) and gi±{z) = 53^1(2) for z G [1, +00); 

• g2 and 54 are analytic on C \ (—00, —1] and g2±{z) = gi^{z) for z G (—00, —1]; 

• as z — >■ 00, 

51 (z) = -i^{-z) - (1 - a/a^)z - \{-z)-^'^ + O (z'^/^) , 
g2{z) = -i,{z) + (1 - a/a^)z - ^^"^2 + o (z-^/^) , 
g^{z) = ipi-z) - (1 - a/a2)z + ^(-^)-^/' + O (z'^/^) , 
54(-2) = i^iz) + (1 - (7/a2)z + + o (2-3/2) _ 

We define M(3) by 

(6.4) M(3)(2;a) = EoM^^\z;a)dia.g (g-^'^i^^), e-"'^^^^), e-^'^^W^ e-a^s*^) ^ 
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where Eo is the constant matrix 



Eq = I4+ ^-E^S,! — ^-^4,2- 



Then M^^^ satisfies the foUowing RH problem with S 



Lemma 6.3. The function M^^^ as defined in (6.4) has the foUowing properties 

(1) M^^^ is analytic for z G C \ Sjy'(3); 

(2) Mf (z) = m1')(z)J,,,(3), for z e 2^,(3)/ 

(3) As z ^ 00 we have 

M^^\z; a) = {l + (z-^)) diag ((-z)-V4^ ^-i/4^ (_^)i/4^ ^i/4^ 

The jump matrices Jj^,/ (3) are given by 



X J^^{2) diag ( e 



And more explicitly, 

•^M(3)(^) 



/I 0\ 

1 

10 

\0 -1 0/ 



, zG (-00,-1), Jm(3)(^) 



/ 1 0\ 

10 

-10 

V 1/ 



, zG(1,oo), 



Jm(3) (^) 



/I \ 

10 e°'^(^2~^4) 

1 



•^M(3)(^) 


















1 






A 







(91- 


-93) 









1 











= 










1 








Vo 












/I 





















1 






















1 







Vo 




(94- 


-92) 











/ 


1 










°\ 











1 












(93- 


^91) 





1 







v 













1/ 






1/ 



ZE (-1,0), 



ZG(0,1), 



, arg(z + 1) = ±(7r - (/7i), 



arg(z - 1) = ±ipi, 
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J 



M{3) 



( 1 

(92-91) 










1 









Y -|-ga^(93-94) 

1 



argz = ±(/92, 



and, finally, 



J 



M(3) 



(^) 



/I 




Vo 



-|-ga^(9l-92) 
1 









1 

a^{94- 



93) 





1/ 



argz = ±(7r - Lp2)- 



Let us discuss the limiting behavior of the jump matrices Jj^iz) as a ^ +00. First assume 
that £7 = 0. In that case, it is not hard to check from (6.3) that in the hmit a — ?■ +00, 
the jump matrices on the contours other than (—00, —1) and (l,oo) are exponentiahy close 
to the identity. To have exponential decay in the upper left block of the jump matrices on 
Ffc, k = 2,3,7,8 we also need 992 ^ vr/G. This can always be achieved deforming these rays 
if necessary. The exponential decay holds uniformly, when we stay away from and ±1, 
suggesting the construction of a global parametrix that solves the limiting RH problem and 
local parametrices near and ±1. Now, if o" 7^ but a large, then near the origin it may 
happen that we have exponential growth of some of the entries of J^,f{3) instead of decay. 
However, this only happens very close to the origin where we are going to construct a local 
parametrix anyway. 

Summarizing, in the next steps we will construct first a global parametrix and then local 
parametrices around and ±1. For our purposes, the local parametrix at the origin is most 



important and it is precisely here that RH problem 2.5 for the Hastings-McLeod solution to 
the Painleve H equation appears. 

6.4 Construction of the global parametrix p(°°) 

In this part we will construct an explicit solution to the RH problem for M'^^^ after taking 
the limit a — )• 00. That is, we want to solve the following RH problem. 
We look for a 4 x 4 matrix- valued function P^^\z) satisfying 

(1) P(°^)(z) is analytic for z € C \ ((-oo,-l] U [l,oo)); 

(2) p(°^) has the following jumps 



(6.5) 



(6.6) 



p|°°)(x) = p!°°)(x) 








1 







1 








-1 











Vo 








1/ 


A 



















1 








1 





Vo 


-1 





0/ 



X G (1, 00), 



X e (—00, — 1) 



where the contour is oriented from left to right; 
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(3) As z — )• oo we have 

p(-)(z) = {l + (z-i)) diag {-z)^/\ z^/^) A, as z ^ oo. 

It is straightforward to check that 

p^^\z) = diag ((1 - (1 + (1 - (i + a 

solves the RH problem. 

6.5 Local parametrices p(='=^) around ±1 

Local parametrices P^^^^ in a neighborhood Ui of 1 and a neighborhood of —1 can be 
constructed using Airy functions. This construction is well-known in the literature and will 
therefore be left to the reader. Also, we will not need the exact form of the solution for our 
analysis. 

6.6 Local parametrix P*^*^^ around 

In this part we will construct a local parametrix around zero out of solutions for two 2x2 
size model RH problems. As can be expected from the structure of the Riemann surface on 
the critical curve r = -\/a+~2, one of these will be related to the PII model RH problem (RH 



problem 2.5), while the other one is constructed from the sine function. More precisely we 



use the following model RH problems. 

Rotated PII model RH problem The first model RH problem will be nothing more 



than a rotated version of RH problem 2.5 It is immediate that 

$(C) = ^{iO, 

satisfies the following RH problem. 

RH problem 6.4. We look for a 2 x 2 matrix-valued function ^{C) satisfying 

(1) $(C) is analytic /or C G C \ S^; 

(2) $+(C) = ^-(C)J5.,/orCGS5,; 

(3) As ( ^ oo we have 



e 3 





(4) ^'(C) is bounded near C = 0. 
The contour consists of the rays arg((^) = ±tt/3 and arg(C) = ±27r/3. The jumps on 



these rays are presented in Figure 15 
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Figure 15: The jump contour in the complex (-plane and the constant jump matrices 
on each of the rays. 




Figure 16: The jump contour Se in the complex (^-plane and the constant jump matrices on 
each of the rays. 



Rotated sine model RH problem The second model RH problem that we need is almost 
a trivial one. It is related to the sine kernel, but since also here we use a rotated version 
exponential functions arise. 

RH problem 6.5. We look for a 2 x 2 matrix-valued function ©(C) satisfying 

(1) 0(C) is analytic for ( £ C \ Se; 

(2) e+(C) = e_(C)J0,/orCGSe; 

(3) As ( ^ oo we have 

e(c) = (/ + o(c-^))('o' °); 

(4) @{C) is bounded near C = 0. 

The contour Se consists of the rays arg(C) = ±7r/3 and arg(C) = =b27r/3 and the jumps Jq 
on these rays are shown in Figure 16. 
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Figure 17: The jump contour and jump matrices for the local parametrix P^^^ (z) around zero. 



On the circle the matching condition (6.7) is imposed. 



This RH problem has the solution 



0(C) = 








, |argC|<^, 



< largCl < f , 



2n 



< I argCI < vr, 



as can be seen by direct verification. 



Structure of the local parametrix We will define the local parametrix P^^\z) on a disk 
Uq of fixed radius e centered at the origin such that it is analytic inside this disk, except on 



the contour shown in Figure 17 On this contour the parametrix needs to have the jumps 
indicated in the same figure. Moreover, on the circle of radius e we impose the matching 
condition 



(6.7) P(0)(z) \z) = I + 0{l/a), 

We look for a solution to this problem of the form 

p{0)(^\ - p(oo)(s f^iafi{z);iy{z)) 



a — >• oo. 



1. / ^3 91-92 „3 91-92 „3 93-94 „3 93-94 \ 

X diag ( e 2 , e 2 , e 2 , e 2 j . 
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We now define the conformal maps /i and /2 and the map u, starting with the first. We 
define 



{I - zf/^ - {1 + z) 



3/2 



1/3 



for z in a small neighborhood Uq of the origin. A simple Taylor expansion shows that the 
cubic root can be taken such that /i is a conformal map satisfying 



Then we define v by 



/i(z) = 2-5/3z + 0(z3), 



as z 



0. 



az 



It is clear that we can choose Uq sufficiently small such that z i— )• is analytic in Uq. Note 
that 1/(0) = 2^/^ a. By the definition of gj in (6.3) and this choice of /i and v, we have 



(6.8) 



a\h{zf - au{z)h{z) = a\g2{z) - gi{z))/2. 



Finally, we define /2 by 



h{z) = \{gA{z) -gz{z)). 



Note that a simple Taylor expansions shows that 



(6.9) 



9z{z) - gdz) = -4z + ^z^ + 0{z'') 



as z 



0. 



Defined in this way, and possibly after some slight deformation of contours in the rotated PII 
and rotated sine RH problems, p('^)(z) already satisfies the jump conditions. It then only 



remains to check the matching condition (6.7) 



Matching condition Let us fix z on the boundary OUq of the disk with radius e. Then 
afi(z) and a? f2{z) both tend to infinity as a — t- c«. We may therefore plug in the asymptotic 



formula for ^ and 0. After a quick calculation using (6.8) it turns out that 



T / /. / \ / \\ 1- ( o3 31-32 3 31-92 

^{ah{z)-v{z))dmg[e " 2 , 2 



e(a3/2(z))diag e 



„3 93 -g4 „3 93-94 
" 2 e 2 



h + 0{l/a), as a — )• +(x>, 
12 + 0(1/0), as a — )• +00. 



It is then clear that also the matching condition (6.7) is satisfied. 
6.7 Final transformation: M^^^ H- M^^) 



Finally we define a last matrix function M^^\z; a) as 

'M(3)(z;a) (PM)"' 
M(3)(z;a) (P(0))"' 



(6.10) 



M(^)(z;a) 



M(3)(z;a) 



for z G C \ {Uo U f/i U U-i U ^M(3)), 
for z eUo\ Sjvf(3), 
for z eUi\ Sjvf(3), 



^M(3)(z; a) (-P^^^)) ^ for z G U-i \ S 



M{3) • 
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This function is analytic in C \ 'Sj^,j(4) , where 5]jyj(4) is obtained from S^^{3) by removing the 
rays (—00, —1] and [1,cxd), adding three smah circles around —1, 0, and 1 and finally deleting 
the parts of 5]^^(3) within these circles. Then all jumps on Sj^,j(4) tend to the identity matrix as 
a —7- +00, both uniformly and in L^-sense. Indeed, for the jump on dUo this is a consequence 



of the matching condition (6.7). Likewise the jumps on dUi and dU-i tend to identity if 



one constructs the right local parametrices around ±1. Moreover M^^^{z; a) is normalized at 
infinity, meaning that 

M^^\z;a) = I + 0{l/z), as z ^ 00. 
Then the following result follows by standard methods. 
Lemma 6.6. As a — )• 00 we have 



M'-^\z;a) =1 + 



1 



a(l + \z\ 



uniformly for z £ C \ Sjv/(4) . 



This is the final goal of the steepest descent analysis and the starting point of the proof 
of the main theorem. 

6.8 Proof of Theorem [2:6] 



Proof of Theorem 2.6. We only give the proof for the case m, f > 0. Other cases can be proved 



similarly. First note that by (|2.15|), (|2.13|), and taking the transpose we obtain 

Kcr{u,V]S,t) 



2TTi{u — v) 



(-1 1 M-^ {iu-s,t)M {iv]s,t) 



/1\ 
1 



Vo/ 



Fix u,v > ^ and let a be a large positive number. Using (6.1) we write 

1 



Kci {a^u, c?v\ ^ — a(l — a jc?^ 



2TTi{u — v) 



-1100) 



X m(^) (ia\; a) ^M^^) {ia^v; a) 
By the transformations M^^) ^ M^^) M^^) in ([6^ and K^ we obtain 



/1\ 

1 





Kcr [a^u, a^v; a, |a^, — a(l — cr/a^)) 



1 



27ri{u — v) 



-110 0) 



X M(3)(i^;;a)diag (^e^'9iM^^a3g2(iv)^^a^g3(iv)^^a3g^iiv)'^ 



/1\ 

1 
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In the next step we unfold the transformation AI^^^ i— t- M^^^ in (6.10). Assume that u,v £ Uq. 
Then we obtain 



^a^(gi{iu)+g2{iu))/2 
pa^{gi{iv)+g2(iv))/2 



Kcr {c?u, a^v; ^a^, — a(l — cr/a^)) 



1 



2-Ki{u — v) 



(-1 1 Q^^h{^h{^n),u{iu))-^ 0V(-)(zu)-iM(4)(^n;a)-i 



X M(^)(w;a)p(°°)(ro) 



■^{afi{iv)-u{iv)) 




/1\ 

1 





Now we make the change of variables u = and v = We consider x and y as 

fixed so that u, v — )• as a — )• +oo. Note that after this change of variables 



afi(iu) — ix, 
afi{iv) iy, 



as a — )• +00, 
as a — )• +00. 



Also note that i^{iu) i^{0) = 



It follows from standard considerations (see also (4.44)) that 



M(^) {iu; ay^M'-^') {iv;a) = I + 



u — V 
a 



I + O 



x-y 
a2 



as a — )■ +00, 



uniformly for x and y in a compact subsets of M. Observe also that 

P^°°\iv) = {I + 0{l/a))A, asa^+oo, 
P^°°\iu)-^ = A-^{I + 0{l/a)), as a ^ +oo. 



Therefore, 



(6.11) 



lim 



a— >+oo 



ga3(gi(m)+g2(m))/2 
^wi{gi{iv)+g2{iv))/2 
1 



Kcr (2^'^ ax, 2^/^ay; \ , — a(l — (T /a) 



1 1) *(ix; 2^'^a)-^^{iy- 2^/^) [J 



2iii{x — y) 

= 2^1(^3^ (-1 i)*(-;#'M-'*(-.;2V3,)(;). 

The convergence is uniform for x and y in compact subsets of M. Now observe that ^{z) 
satisfies the symmetry 

\ J) ^(.; 2^3^) (J J)=^(-.;2V3,). 

This can be easily verified by checking that both sides solve the same RH problem and noting 
that the RH problem has a unique solution. Finally, after inserting this symmetry in (6.11) 
and using (2.19) we obtain 

^ci?(g'i_(iu)+g2{iu))/2 



hm 2'>l\'' 



a^{gi(iv)+g2{iv))/2 



Ker (2^/3„x, 25/3ay; , -a(l - a /a)) = Kpu{x, y; 2^/^). 



This completes the proof of Theorem 2.6 for u,v > 0. 



□ 
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A Asymptotic behavior of Kcr and ii'tac along the diagonal 



Here we prove the asymptotic expansions (2.21) and (2.22) for the diagonal of Kcr as defined 
in (2.15) and Ktac in (2.20) . First we note that that the diagonal terms can be obtained by 
taking the Umit u — )• v in (2.15) and (2.20), which leads to 

(A.l) 



(A.2) 



Kcr{u,u) 



Kta.c{u,u) 



— (1-10 Q) M-\tu) — {iu) 



^(-10 1 0)M-(.)^(.) 



1 


Vo/ 



1 



for u > 0. Next we write 



M(c) = nc)Q(c), 

where Y has a full and uniform asymptotic expansion in inverse powers of C, 



(A.3) 
and 



y(C) = h + ^ + [C^) , as c ^ oo, c G c, 



Q(c) = diag((-c)-l/^ r (-c)'/^ c^/^) 



Here V'l and '(/'2 are as defined in (^2.9|-(2.10). A straightforward computation shows 



Q'(C) = i?(C)Q(C), 



where 



( + ' 



\ 



«r2 



IS2 



\ 





-iriC — is 



4C 



+ t 



4C 



/ 



Then 



dM 



M-\0—{0 = Q-\c) {Y-HOY'iO + ^(0) Q(C). 



When we let C tend to oo the contribution of the first term is of order 0{C ^) so we will 
concentrate on the second term. For the tacnode one-point correlation function we find 

(A.4) 

n\ 

Kts,ciu,u;ri,r2,si,S2) = — {-I 1 O) Q^^{u)R{u)Q. 



2m 



1 




1 

Vo/ 



27ri V 4n V 



TT TT-v/u 47rn 
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as u +00. Note that there is no constant term in the expansion and that the 1/u has a 
highly oscihatory coefficient of modulus 1. For the two-matrix model we can compute 



Kcr{u,u;s,t) = — (1 -1 0)Q~'^{iu)R{iu)Q{iu) 



(A.5) 



/1\ 
1 





( (_i^)l/2 + (i^)l/2 + 2i + s(-i«)-V2 + s(m)-V2) + 0(^-3/2) 



t 

+ - + 



as It —7- +00. Here we do have a constant contribution, but there is no 1/u term. Therefore, 
these two expansions cannot be the same and hence the kernels must be different. 
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